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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

EXCLUSIVE  MEASUREMENTS  IN  5 D*NNX 

By 

Antonio  I.  Rubiera 
August  2000 

Chairman:  J.  Yelton 
Major  Department:  Physics 

We  report  the  first  observation  of  exclusive  decays  of  the  type  B — > D*  NNX , 
where  N is  & nucleon.  Using  a sample  of  9.7  x 10®  BB  pairs  collected  with  the 
CLEO  detector  operating  at  the  Cornell  Electron  Storage  Ring,  we  measure  the 
branching  fractions  B (B^  — >■  D*~pp7r'^}  = ±1-0)  x 10“^,  and  B(B^-^ 

D*~pn  ) = (14.5l|;o  ± 2.7)  x 10““*.  The  charge  conjugate  process  is  implied  in 
the  reconstruction  of  B^—^  D*~  p p tt"*".  However,  in  the  reconstruction  of  R®— > 
D*~  p n,  only  the  mode  with  the  antineutron  is  used  in  our  measurement  because 
neutrons  do  not  have  a distinctive  annihilation  signature. 

Antineutrons  are  identified  by  their  annihilation  in  the  Csl  electromagnetic 
calorimeter.  Since  we  are  unable  to  isolate  a sample  of  antineutrons  in  data,  we 
use  antiproton  annihilation  showers  in  a A — > sample  to  define  the  antineu- 

tron selection  criteria.  We  find  a discrepancy  for  antiproton  annihilation  showers 
between  the  Monte  Carlo  and  data,  which  we  assume  affects  antineutrons  as  well. 
We  increase  the  raw  yield  for  R®— >■  D*~  p n by  21%  to  correct  for  this  discrepancy. 

The  possible  contributions  from  J5®  D*~  Df  with  Df  —^pn  and  5®  D* 

with  £)*+  Dfj  and  Df  ^ pn  are  eliminated  from  the  analysis  by  rejecting 

• • 
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events  with  1.91  GeV  < Mp+fi  < 2.04  GeV  for  a loss  of  9%  in  the  reconstruction 
efficiency.  We  fail  to  find  evidence  for  the  decay  Df  —^pfi. 

We  search  for  possible  contributions  to  the  resonant  substructure  of  )■  D*~ 
p n and  D*~  p p 'k'^  due  to  a heavy  charmed  baryon  decaying  strongly  to  p 
D*~  for  D*~  p p TT'^  and  n D*~  for  D*~  p n,  as  well  as  a resonance  of 
the  virtual  W decaying  to  pp'ir'^.  We  also  study  the  possible  effect  of  feed-down  A 
baryon  contributions  to  the  background  for  both  modes,  as  well  as  the  5^— > D*~ 
p p TT'^  signal.  No  conclusive  evidence  is  found  for  a measurable  contribution  from 
the  aforementioned  contributions  to  the  resonant  substructure. 

Antineutrons  are  used  for  the  first  time  in  the  exclusive  reconstruction  of  a 
B meson.  By  finding  conclusive  evidence  for  the  existence  of  decay  modes  of  the 
type  B DNNX,  we  challenge  the  assumption  that  the  B — > Baryons  rate  is 
dominated  by  decays  to  charmed  baryons. 


Vlll 


List  of  Figures 


1.1  A Feynman  diagram  for  neutron  beta  decay  in  Fermi  Weak  Theory 

1.2  A Feynman  diagram  for  neutron  beta  decay  in  the  Standard  Model 

1.3  A Feynman  diagram  for  D*  D^-Ksoft 

1.4  A second  Feynman  diagram  for  D*  — > D^iTsoft 

1.5  A color-allowed  Feynman  diagram  for  B~  D°7r“  for  one  quark 

color 

1.6  A color-suppressed  Feynman  diagram  for  B — > for  one  quark 

color 

1.7  A Feynman  diagram  for  B~  — )■  Ac^pn~ 

1.8  Two  Feynman  diagrams  for  D*~  p n 

1.9  A Feynman  diagram  for  B^^  D*~  p p tt"*"  


9 

9 

11 

12 

14 

15 
22 
26 
27 


2.1  Cross  section  into  hadrons  from  the  collision  of  e"''e  beams  at  CESR 
as  measured  by  the  CLEO  II  detector  in  the  energy  range  9.44  GeV 

to  10.62  GeV 

2.2  Front  view  of  the  CLEO  II  detector  

2.3  Side  view  of  the  CLEO  II  detector  

2.4  dE/dx  vs.  track  momentum 

2.5  Main  drift  chamber  (DR)  wire  arrangement  

2.6  Time-of-Flight  vs.  track  momentum  

2.7  Layout  of  CLEO  II  detector  showing  barrel  and  endcap  calorimeter 

sections  


29 

31 

32 

38 

39 

41 

42 


3.1  m£)*-m£)0  in  GeV  for  B D*X  with  in  CLEO  II  ... 

3.2  m£)*-W£)0  in  GeV  for  B — > D*X  with  in  CLEO  II. 5 . . 

3.3  mo*-m£)0  in  GeV  for  B D*X  with  tt®  in  CLEO  II  . . 

3.4  rriD*-'^D°  ™ GeV  for  B — )■  D*X  with  tx^  in  CLEO  II. 5 . 

3.5  Tnjr)*-m£)o  in  GeV  for  B — )■  D*X  with  D^—^K^tx  tx'^tx  in  CLEO  II 

3.6  mD*-mDQ  in  GeV  for  B D*X  with  D^—^K'^tx  tx'^tx  in  CLEO  II. 5 

3.7  Emain  VS.  PQCD  for  protons  and  antiprotons  in  CLEO  II  

3.8  Emain  VS.  PQCD  for  protons  and  antiprotons  in  CLEO  II  

3.9  Emain  In  GeV  without  proton  requirement  for  n’s,  y’s  from  tt^’s,  and 

KEs . . . . 

3.10  Emain  ln  GeV  with  proton  requirement  for  n’s,  y’s  from  7t°’s,  and 




57 

58 

59 

60 
61 
62 

70 

71 

74 

74 


IX 


4.1  AE  vs  Mbc  distribution  for  > D*~  p p in  CLEO  II  ON 

resonance  data 

4.2  AE  vs  Mbc  distribution  for  J5°— > D*~  p p in  CLEO  II. 5 ON 

resonance  data  SI 

4.3  Mbc  (in  GeV)  for  D*~  p p -k+  in  CLEO  II 82 

4.4  Mbc  (in  GeV)  for  D*~  p p in  CLEO  II.5 83 

4.5  Mbc  for  D*~  p p in  data 86 

4.6  AE  (in  GeV)  for  > D*~  p p -k^  in  CLEO  II/II.5  data  87 

4.7  7Ts  from  D*  momentum  for  B^-^  D*~  p p in  data  and  Monte 

Carlo  S8 

4.8  Mbc  in  data  and  generic  Monte  Carlo  for  S'* pp  (in  GeV)  92 

5.1  A Feynman  diagram  for  Df  p 96 

5.2  Mp+n  in  GeV.  B^  Df  D*~  with  Dj^pn 97 

5.3  Mp+fi  (in  GeV)  from  a reconstruction  of  5°->  D*~  p n in  a signal 

Monte  Carlo  B^  -4  D*+  D*~  with  Df  p n 98 

5.4  (white)  is  5°  Df  D*~  with  p n.  (dashed)  is 

D*~  with  Df  p n.  Mp+n  (in  GeV)  99 

5.5  (white)  is  B^  — >■  D'^  D*~  with  Df  p n.  (dashed)  is  B^  -4 

D*~  with  p h.  mBO  (in  GeV)  99 

5.6  (white)  is  inclusive  of  all  contributions,  (dashed)  is  after  exclusion 

of  both  -4  p n contributions  191 

5.7  (white)  is  ON  resonance,  (solid)  is  OFF  resonance  102 

5.8  Mbo  for  B‘^-4  D*~  p n in  CLEO  II  ON  resonance  data  103 

5.9  Mbo  for  jB“-4  D*~  p n in  CLEO  II. 5 ON  resonance  data  103 

5.10  Mp+n  for  B°-4  D*~  p n (in  GeV)  in  CLEO  II/II.5 104 

5.11  Mp+n  vs.  tubo  (both  in  GeV)  105 

5.12  Inclusive  tubo  in  data  for  B°-4  D*~  p n (in  GeV)  106 

5.13  niBO  in  data  for  B°-4  D*~  p h in  GeV  excluding  both  D+  ->  p n 

contributions  ' 

5.14  Antineutron  momentum  distribution  (in  GeV) 112 

5.15  tubo  in  GeV  for  D*~  p n selection  criteria  applied  to  selected 

signal  Monte  Carlo  background  modes  114 

5.16  tubo  for  B^^  D*~  p n in  GeV.  (white)  is  data,  (solid)  is  generic  MC  115 

5.17  (in  GeV)  for  neutrons  and  antineutrons  in  B°-4  D*~  p n 

signal  Monte  Carlo 


X 


List  of  Tables 


1.1  Quark  families  

1.2  Anti-quark  families 

1.3  Estimated  quark  masses 

1.4  Selected  meson  masses 

1.5  Selected  baryon  masses  

1.6  Lepton  masses  

1.7  Some  characteristics  of  bosons  

1.8  Some  B meson  decay  mechanisms 

1.9  Representative  B — >■  Baryons  decay  mechanisms 

2.1  Some  e'^e~  annihilation  final  states 

2.2  Vertex  detector  and  drift  chamber  parameters 

2.3  Drift  chamber  resolution 

2.4  Momentum  resolution  {Spt)  for  the  CLEO  II  data  at  selected  values 

of  

2.5  Radiation  lengths  of  material  before  barrel  section  of  the  calorimeter 

2.6  Energy  and  angle  resolutions  for  a photon  in  the  barrel  at  two  values 

of  cluster  energy  

3.1  The  CLEO  data  

3.2  Particle  identification  criteria 

3.3  mpo  and  mu^-rriDO  double  Gaussian  data  cuts 

3.4  D*  yield:  

3.5  D*  yield:  

3.6  D*  yield:  D°-^R+7r“7r'^7r~  

3.7  Shower  types  and  energy  measured  in  calorimeter 

3.8  Antineutron  shower  selection  criteria 

3.9  Shower  population  in  a generic  BB  Monte  Carlo  sample  without 

proton  requirement ^ 

3.10  Shower  population  in  a generic  BB  Monte  Carlo  sample  with  proton 

requirement 

4.1  AE  fit  results  for  B°— > D*~  p p signal  MC  in  CLEO  II  . . . . 

4.2  AE  fit  results  for  BM  D*~  p p 7t+  signal  MC  in  CLEO  IL5  . . . 

4.3  Mbc  distribution  fit  results  for  B®— > D*  p p signal  MC  in 

CLEO  II  


4.4  Mbc  distribution  fit  results  for  D*  p p signal  MC  in 

CLEO  II.5  79 

4.5  Results  in  the  CLEO  II/II.5  data  of  various  fits  to  Mbc  for  B^—^ 

D*~  p p TT'^  84 

4.6  B^—^  D*~  p p : number  of  events  found  per  mode 85 

4.7  % composition  of  generated  (phase  space  MC)  and  events  recon- 
structed (data)  of  Mp^+  and  Mp^+  90 

4.8  Estimate  of  systematic  uncertainties  (in  %)  for  B^-^  D*~  p p tt"*"  . 93 


5.1  Antineutron  shower  selection  criteria 

5.2  B^-^  D*~  p n signal  MC  e^c’s  and  widths  in  CLEO  II 

5.3  B^^  D*~  p n signal  MC  cmc’s  and  widths  in  CLEO  II.5  

5.4  Measurements  of  B^—^  D*~  p n branching  fractions  

5.5  Results  of  various  fits  to  m^o  for  B^—¥  D*~  p n with  or  without 

the  Dj"  p h contributions 

5.6  A selection  criteria 

5.7  A double  Gaussian  fixed  parameters  _ 

5.8  eannihiiation-  (%)  of  match  1 or  2 antiprotons  in  A’s  passing  annihi- 
lation shower  selection  cuts  for  six  momentum  spectra  

5.9  Correction  Factor  (C.F.)  as  a function  of  Pproton  and  % of  antineu- 

trons found  in  signal  Monte  Carlo  for  each  momentum  range  in 
CLEOII/II.5 

5.10  Neutron  shower  selection  criteria  

5.11  Estimate  of  systematic  uncertainties  (in  %)  for  B^—¥  D*  p n . . . 


94 

100 

100 

107 

108 
no 
no 

111 


113 

116 

118 


6.1  Exclusive  measurements  in  R — > Baryons  to  date 


119 


Xll 


CHAPTER  1 
INTRODUCTION 


The  study  of  the  elementary  particles  at  a wide  range  of  interaction  energies 
has  occupied  scientists  since  the  discovery  of  the  electron  in  1896  [1].  Particle 
physics  has  evolved  from  a field  involved  in  the  discovery  of  new  particles  to  one 
devoted  to  their  systematic  study.  A logical  structure,  currently  explained  by  the 
Standard  Model  of  elementary  particles  [2,  3,  4],  has  been  unveiled. 

The  Standard  Model,  however,  offers  an  incomplete  picture  of  some  experi- 
mental high  energy  physics  results.  The  results  we  describe  are  amongst  these. 
Despite  some  weaknesses,  the  Standard  Model  has  withstood  intense  experimental 
scrutiny,  and  while  some  of  the  equations  are  currently  not  calculable,  evidence 
has  not  been  found  for  physics  beyond  the  Standard  Model.  Particles  acquire  their 
masses  in  the  Standard  Model  via  the  Higgs  mechanism,  that  requires  the  existence 
of  a massive  gauge  boson,  the  Higgs  boson,  that  has  yet  to  be  found. 

We  first  introduce  particle  physics  phenomenology  -particles  and  their  decays, 
which  is  the  information  most  often  used  by  the  practitioner  of  experimental  high 
energy  physics.  A discussion  of  the  current  status  of  B physics  phenomenology 
follows.  Finally,  we  review  B Baryons  previous  to  our  results,  concluding  with 
an  overview  of  our  work  and  its  significance. 
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1.1  Matter 

Matter  is  composed  of  three  types  of  indivisible  constituents;  leptons,  gauge 
bosons,  and  quarks.  Leptons  and  gauge  bosons  are  found  individually  in  Na- 
ture. Quarks  combine  in  two  types  of  arrangements  to  form  hadrons.  The  first 
arrangement  is  of  the  form  quark-antiquark,  and  is  called  a meson.  The  second 
arrangement,  three  quarks,  is  called  a baryon.  Mesons  and  baryons,  collectively 
known  as  hadrons,  comprise  all  the  known  forms  of  matter  consisting  of  quarks. 

Hadronic  matter  is  said  to  be  colorless.  A color  is  assigned  to  each  quark  in  a 
meson  or  a baryon.  The  three  quarks  in  a baryon  each  have  a different  color— red, 
green,  blue—  and  the  combination  of  all  three  colors  yields  a colorless  hadron.  The 
quark  and  the  antiquark  in  a meson  form  a color-anticolor  pair  (e.g.  QbiueQbl^)- 

1.1.1  Hadrons 

There  are  three  families  of  quarks,  each  consisting  of  two  types  of  quark:  an 
up-type  quark,  with  electromagnetic  charges  -t-2/3  that  of  the  electron,  or  (-1-2/3) 
9e,  and  a down-type  quark,  with  (-1/3)  Qe-  Every  type  of  quark  is  called  a flavor 

of  quark. 

The  first  and  lightest  family  consists  of  the  up  {u)  and  down  (d)  quarks.  The 
next  family,  with  heavier  quarks,  is  the  charm  (c)  and  strange  (s)  family.  Even 
heavier  still  is  the  third  family:  the  top  (t)  and  bottom  (b)  quarks.  All  are  shown 

in  Table  1.1. 

For  every  matter  constituent  there  is  an  anti-matter  constituent  with  opposite 
electromagnetic  charge  and  equal  mass,  as  shown  in  Table  1.2. 
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Table  1.1:  Quark  families 


up-type  (+2/3)  Qe 

up:  u 

charm:  c 

top:  t 

down-type  (-1/3)  qe 

down:  d 

strange:  s 

bottom:  b 

Table  1.2:  Anti-quark  families 


up-type  (-2/3)  Qe 

u 

c 

t 

down-type  (+1/3)  Qe 

d 

s 

b 

Quarks  are  not  found  individually  in  Nature.  Their  masses  can  be  estimated 
by  the  mass  spectrum  of  mesons  and  baryons  that  has  been  measured  to  date. 
In  Table  1.3  we  show  current  best  estimates  of  the  lower  and  upper  limits  for 
quark  masses.  The  estimated  masses  for  quarks  are  useful  in  current  Standard 
Model  calculations.  However,  the  quark  mass  estimates  have  large  uncertainties, 
especially  in  the  case  of  the  up  and  down  quarks.  The  values  quoted  in  this  Table 
and  the  next  values  quoted  in  this  section  are  from  the  1998  Review  of  particle 
physics  [5].  The  mass  of  the  top  quark  mt  = 173.8±5.2  GeV. 


Table  1.3:  Estimated  quark  masses 


Quark 

Low  mass 

High  mass 

u 

1.5  MeV 

5 MeV 

d 

3 MeV 

9 MeV 

s 

60  MeV 

170  MeV 

c 

1.1  GeV 

1.4  GeV 

b 

4.1  GeV 

4.4  GeV 

The  proton,  for  example,  has  the  quark  content  (uud),  and  one  of  the  pions, 
the  TT'^,  has  the  quark  content  {ud).  Each  type  of  quark  is  considered  a distinct 
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flavor.  The  group  of  mesons  containing  a charm  quark  and  a light  antiquark  (one 
of  s,  u,  d)  is  called  the  charm  mesons. 

Meson  and  baryon  masses  are  known  to  varying  degrees  of  accuracy,  as  shown 
in  Table  1.4  for  selected  mesons,  and  Table  1.5  for  selected  baryons. 


Table  1.4:  Selected  meson  masses 


Meson 

Quark  content 

Mass 

7T^ 

ud 

139.56995  ± 0.00035  MeV 

su 

493.677  ± 0.016  MeV 

cu 

1864.6  ± 0.5  MeV 

bd 

5279.2  ± 1.8  MeV 

Table  1.5:  Selected  baryon  masses 


Meson 

Quark  content 

Mass 

proton 

uud 

938.27231  ± 0.00028  MeV 

A 

sud 

1115.683  ± 0.006  MeV 

Ac 

cud 

2284.9  ± 0.6  MeV 

Afe 

bud 

5624  ± 9 MeV 

1.1.2  Leptons. 

The  leptons,  the  e,  the  /i,  and  the  r,  are  fundamental  particles.  Each  can  be 
produced  with  an  accompanying  neutrino,  z//,  with  / = e,  /i,  or  r.  Neutrino  physics 
in  the  near  future  will  yield  a better  understanding  than  is  currently  available. 
Neutrinos  have  been  thought  to  be  massless  and  not  to  mix  (that  is,  each  neutrino 
flavor  was  thought  only  to  interact  with  its  lepton  partner).  The  recent  observation 
of  neutrino  mixing,  however,  implies  that  neutrinos  have  mass  [6].  Unlike  quarks. 
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leptons  are  observed  as  single  particles  with  well-measured  masses.  In  Table  1.6 
we  show  the  current  mass  measurements.  Large  differences  in  masses  are  found 
between  the  electron,  the  muon,  and  the  tau. 


Table  1.6:  Lepton  masses 


Lepton 

Mass 

electron: 
muon;  jj,'^ 
tau: 

0.51099907  ± 0.00000015  MeV 
105.658389  ± 0.000034  MeV 
1777.05l^;i  GeV 

Helicity  is  the  orientation  of  a particle’s  momentum  vector  in  respect  to  its 
spin.  Helicity  can  be  right-handed  or  left-handed  for  all  fermions  except  neutrinos, 
which  are  only  left-handed.  Antineutrinos,  likewise,  are  only  right-handed. 


1.1.3  Gauge  Bosons 

The  mediators  of  the  physical  forces,  the  particles  that  allow  decays  to  take 
place,  are  called  gauge  bosons.  In  Table  1.7  we  outline  their  masses  and  the 
types  of  interaction  that  each  mediates.  I = e,ii,T,  n=up-type  quark(u,  c,  t),  and 
d=down-type  quark(d,  s,  b). 


Table  1.7;  Some  characteristics  of  bosons 


Boson 

Force 

Mass 

Decays 

photon 

gluons 

electromagnetic 

strong 

weak 

weak 

massless 

massless 

80.41  ± 0.10  GeV 
91.187  ±0.007  GeV 

7 ^ l+l~ 

9 ^ QQ 

■W+  l+vi,ud 

W+  -4  in-,vu,q^q~ 
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1.1.4  Spin  and  Statistics 

Matter  is  also  characterized  by  the  statistics  obeyed.  Leptons  and  quarks  are 
fermions,  obeying  Fermi  statistics,  an  example  of  which  is  the  Pauli  exclusion 
principle  for  electrons  occupying  the  same  shell  in  an  atom.  The  spin  of  leptons 
and  quarks  is  (±|).  Bosons  obey  Bose  statistics,  which  allow  an  infinite  number  of 
particles  to  occupy  the  same  energy  state.  Bosons  have  integral  spin  (0,  ±1).  The 
diflferent  relative  alignments  of  the  spins  of  the  individual  quarks,  together  with 
the  addition  of  angular  momentum,  results  in  a large  number  of  possible  states. 


1.1.5  The  CKM  Matrix 

The  interactions  of  quarks  from  different  flavor  families  are  suppressed  with 
respect  to  those  within  the  same  family.  In  order  to  correct  for  this  discrepancy, 
the  characteristics  of  a given  decay  as  prescribed  by  weak  theory  are  adjusted  using 
the  flavor  mixing  3x3  unitary  and  complex  matrix  V , the  Cabibbo-Kobayashi- 
Mciskawa  matrix  [5],  which  transforms  the  weak  d,  s,  b quark  states  to  the  physically 
measured  d! , s'  ,b'  quark  states,  to  arrive  at  a correct  theoretical  understanding  of 
a weak  decay  for  all  quark  flavors: 


VcKM  — 


^ Kd(0.9745  - 0.9760) 
Kd(0.217  - 0.224) 
I/td(0.004  - 0.013) 


K,(0.217  - 0.224) 
14,(0.9737  - 0.9753) 
14,(0.035  - 0.042) 


I4b(0.004  - 0.013)  ' 

44fc(0.036  - 0.042) 
I4fc(0.9991  - 0.9994)  ^ 


with 
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The  diagonal  elements  of  F are  « 1,  implying  that  decays  which  involve  these 
CKM  matrix  elements  are  Cabibbo-allowed,  whereas  all  other  decays,  which  involve 
off  diagonal  elements,  are  Cabibbo-suppressed.  For  example,  a,  b cVF  , W — > 
cs  decay  has  a much  larger  decay  rate  than  a 6 cW  ,W  cd  decay. 

1.1.6  Symmetries 

A particle  and  its  antiparticle  are  said  to  be  symmetric  under  CPT  transfor- 
mations, where  C is  charge  conjugation,  P is  parity,  and  T is  time  reversal.  The 
inversion  of  parity  acts  like  a mirror  in  the  inversion  of  space  coordinates.  CPT 

symmetry  is  valid  for  all  forces. 

A subset  of  this  symmetry  is  CP.  The  product  of  a charge  conjugation  and  a 
parity  inversion  affect  a particle  the  same  as  its  antiparticle.  When  CP  is  violated, 
a particle  will  prefer  a different  subset  of  decays  than  its  antiparticle.  Such  is  the 
case  for  neutral  kaons.The  and  states  in  weak  theory  are  different  from 
the  physically  observed  strong  states,  the  As,  or  A-short  and  the  Al,  or  A-long, 
each  with  its  respective  antiparticle.  Short  and  long  refer  to  the  respective  decay 
lifetimes.  Charge  symmetry  is  obeyed  implicitly  in  neutral  kaon  decay,  but  parity 
is  violated.  CP  violation  is  evident  in  the  difference  in  decay  rates  of  As’s  and 
Al’s  to  final  states  of  two  and  three  pious,  the  former  being  even  under  parity. 
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and  the  latter  odd  under  parity  [7].  CP  violation  in  B decays  and  will  be  the  focus 
of  many  studies  in  the  near  future. 

1.2  Decays 

All  hadrons,  except  the  proton  and  electron,  are  unstable  and  decay  to  lighter 
hadrons,  leptons,  and  bosons.  The  top  quark  decays  well  before  the  time  required 
to  form  hadrons  (baryons  or  mesons).  Therefore,  the  bottom  group  of  baryons 
(e.g.  At,  with  quark  content  bud)  and  mesons  (e.g.  Be,  with  quark  content  be)  is 
the  group  of  hadrons  with  the  largest  masses  currently  found  in  Nature. 

1.2.1  Weak  Decays 

The  b quark  in  the  B meson  decays  via  the  weak  interaction.  The  b quark 
decay  is  often  accompanied  by  the  strong  decays  of  soft  gluons  which  allow  for  the 
hadronization  of  a large  number  of  final  states.  The  Fermi  Electroweak  theory, 
which  aimed  to  explain  neutron  13  decay,  n as  shown  by  the  Feynman 

diagram  Figure  1.1,  introduced  the  neutrino  to  particle  physics,  serving  as  the 
precursor  of  the  Standard  Model.  A full  theoretical  understanding  of  this  decay 
was  accomplished  by  the  Standard  Model,  in  which  this  decay  is  mediated  by 
the  W vector  boson,  not  present  in  the  Fermi  Electroweak  theory.  The  decay 
^ pe~i7e  is  correctly  described  as  the  quark  level  process  d ^ W u,  followed  by 
W~  e'barue,  as  shown  by  the  Feynman  diagram  in  Figure  1.2. 
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Figure  1.1;  A Feynman  diagram  for  neutron  beta  decay  in  Fermi  Weak  Theory 


Figure  1.2:  A Feynman  diagram  for  neutron  beta  decay  in  the  Standard  Model 
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Examples  of  weak  semi-leptonic  decays  of  mesons,  analogous  to  neutron  (5 
decay,  are  ^ 7r°e+r/e,  ^ K~e^v,,  and  D~e^u,.  The  weak  semi- 

leptonic  decays  of  baryons  follow  a similar  pattern.  These  decays  are  referred  to  as 
semi-leptonic  decays  since  the  final  products  are  a combination  of  a weak  leptonic 
decay  and  quark  hadronization.  Leptonic  weak  decays,  in  which  the  unstable 
particle  annihilates  into  an  l+ui  pair,  examples  of  which  are  7t+  and 

}{+  are  the  most  theoretically  well  understood  type  of  decays,  as  it  lacks 

any  final  state  hadronization. 

1.2.2  Strong  Decays 

The  weak  decays  of  the  heavy  (charm  and  bottom)  quarks  in  mesons  and 
baryons  often  involve  the  secondary  emission  of  strongly-interacting  soft  gluons. 
While  we  understand  the  weak  component  of  these  decays  using  the  Standard 

Model,  the  strong  component  is  not  calculable. 

The  strong  interaction  hadronization  process  is  not  well  understood  theoreti- 
cally when  soft,  or  low  momentum,  gluons  mediate  the  decay.  The  Standard  Model 
is  based  on  perturbation  expansions  which  rely  on  convergence.  A process  involv- 
ing soft  gluons  yields  equations  that  are  no  longer  perturbatively  convergent.  This 
stumbling  block  has  prevented  us  from  understanding  many  details  of  unstable 
particle  decay,  particularly  for  the  case  of  heavy  hadrons,  m which  the  large  mass 
available  in  the  decay  implies  a very  large  number  of  possible  final  decay  products 

from  an  equally  large  number  of  soft  gluons. 

The  D*  meson,  for  example,  which  we  reconstruct  in  this  analysis,  is  a spin 

1 meson  that  decays  via  the  strong  interaction.  Two  possible  Feynman  diagrams 
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for  this  decay  are  shown  in  Figures  1.3  and  1.4.  Our  ignorance  about  soft  gluon 
strong  interactions  forbids  us  from  knowing  the  proportion  of  the  total  decay  rate 
due  to  any  one  diagram. 


Figure  1.3:  A Feynman  diagram  for  D*  — )•  D'^'Ksoft 


Figure  1.4:  A second  Feynman  diagram  for  D*  D^T^soft 
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1.3  B Meson  Decays 

The  B meson,  which  is  the  topic  of  our  study,  can  decay  to  a large  number  of 
lighter  particles  by  various  mechanisms,  some  of  which  are  detailed  in  Table  1.8. 
The  b cW~  type  decays  account  for  more  than  95%  of  the  B decays  that  are 
possible  in  the  Standard  Model.  The  combined  semileptonic  decay  rate  for  the 
three  lepton  families  is  25%,  with  hadronic  decays  accounting  for  almost  all  of 
the  remaining  decay  rate. 


Table  1.8:  Some  B meson  decay  mechanisms 


Quark-level  mechanism 

Sample  decay  mode 

Semileptonic  decay: 
b — y cW  , W — ^ 1 Pi 
Hadronic  decays: 
b cW~ ^W~  cs 
b ->  cW~,W~  — ^ ud 

B 

B DDs 
B -4  Dn 

The  W boson  is  colorless.  Its  decay  to  quarks  constrains  the  color  of  both 
quarks  to  cancel.  The  number  of  Feynman  diagrams  for  a given  B decay  by  the 
color  of  the  quarks.  Whereas  decays  mediated  by  an  external  W boson  allow  for 
any  of  the  three  possible  colors  (for  example,  the  tt"  in  B~  ->  D°7r  , as  shown  in 
the  Feynman  diagram  in  Figure  1.5).  Those  in  which  the  W boson  decays  internally 
limit  the  color  of  all  the  quarks  to  be  the  same  as  that  of  the  parent  meson,  as 
shown  in  the  Feynman  diagram  in  Figure  1.6.  The  latter  quality  is  referred  to  as 
color  suppression.  Decays  that  only  have  an  internal  W boson-mediated  diagram, 
such  as  B^  are  color-suppressed  decays. 
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Figure  1.5;  A color-allowed  Feynman  diagram  for  B — > D^tt  for  one  quark  color 
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Figure  1.6:  A color-suppressed  Feynman  diagram  for  B D^tc  for  one  quark 
color 
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The  decays  b -4-  {c,u)W~  are  mediated  by  the  W~  vector  boson,  while  decays 
of  type  b -4  (s,  d)7  and  b —4  (s,  d)g  are  mediated  by  neutral  bosons.  All  decays 
except  b —4  cW~  contribute  a small  fraction  of  the  total  decay  rate.  However,  it 
is  for  many  of  these  rare  decays  that  our  current  phenomenological  understanding 
allows  for  the  closest  agreement  between  theory  and  experiment. 

It  is  theoretically  allowed,  and  has  been  experimentally  measured,  that  a B 
meson  can  oscillate  to  a before  decay,  allowing  for  neutral  T(4S')  —4  BB  events 
with  either  2 jB^’s,  or  2 B^'s  [8]. 


1.3.1  Quantum  Chromodynamics 


The  Standard  Model  unifies  the  electromagnetic  and  the  weak  interactions.  It 
also  encompasses  the  strong  interactions  in  the  form  of  Quantum  Chromodynam- 
ics (QCD)  [9,  10].  The  strong  coupling  constant  is  smallest  for  short  range 
interactions,  or  for  large  momentum  transfer,  a quality  of  QCD  called  asymptotic 
freedom. 

The  QCD  Lagrangian  is  given  by  [11]; 


\ 

Qred 


^QCD  — i^Qred')  Qblue'i  Q green)  ^ fJi 
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+ h.c. 
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where  the  covariant  derivative  is; 


Dfjt  — d^  + ^igX''  Afi 


17 


|A^  are  the  SU(3)  flavor  matrices,  and  are  the  eight  color  gauge  fields. 
is  the  gluon  field-strength  tensor. 

1.3.2  Heavy  Quark  Effective  Theory 

The  properties  of  the  D and  B mesons,  which  are  composed  of  a heavy  quark 
and  a light  anti-quark,  have  been  described  using  heavy-quark  symmetry  [12,  13] 
by  Heavy  Quark  Effective  Theory  (HQET).  In  a meson,  a quark  and  an  antiquark 
are  confined  to  a bound  state  in  a cloud  of  virtual  quarks  and  gluons  which  need 
to  be  incorporated  into  any  calculation.  In  the  case  of  a heavy  meson,  such  as  a H 
meson,  the  heavy  quark  has  a substantially  larger  mass  than  the  light  antiquark.  In 
HQET,  the  heavy  quark  is  independent  of  the  light  anti-quark.  HQET  assumptions 
simplify  the  Standard  Model  equations,  allowing,  for  instance,  the  comparison  with 
theory  of  experimental  values  of  Vub  and  Vcb- 


The  effective  Lagrangian  that  is  used  to  characterize  B decays  is  given  by 
[14,  15]: 


where  Gp,  the  Fermi  constant,  is  1.17  GeV  and  Jqq  is  the  charged  weak 
current  given  by: 


\ 


\ 


jJ.L  + {uL,CL,tL)'f^VcKM  Sl 


\ 


Jcc  = (^e,  I'n, 
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The  assumption  that  the  heavy  quark  mass  ttiq  is  effectively  infinite  is  used  to 
simplify  the  QCD  Lagrangian.  The  heavy  quark  and  the  light  quark  are  decoupled, 
and  the  effect  of  the  cloud  of  virtual  light  quarks,  light  antiquarks,  and  gluons  is 
assumed  to  be  small  enough  to  be  ignored.  The  QCD  Lagrangian,  of  which  £g// 
is  a simplified  version,  is  simplified  to; 

- mQ)Q 

where  is  the  QCD  covariant  derivative. 

In  the  limit  mg  oo  the  strong  interactions  of  the  heavy  quark  become 
independent  of  its  mass  and  spin  and  the  effective  Lagrangian  is  further  simplified 
to: 

C^HQET  — hyiv  • Dh„ 

where  is  the  effective  heavy  quark  field  and  v is  the  hadron’s  velocity,  which 
is  close  to  that  of  the  heavy  quark. 

In  the  calculation  of  HQET  quantities,  the  strong  interaction  effects  that  are 
non-perturbative  are  grouped  into  a form  factor  that  includes  a dimensionless 
probability  function,  the  Isgur-Wise  function  ^{v  • v)  [16],  where  v and  v are 
respectively  the  initial  and  final  velocities  in  a scattering  or  decay  process.  An 
example  of  the  role  this  function  plays  is  the  elastic  scattering  of  a B meson.  The 
hadronic  matrix  element  for  this  process  is: 

;^{B{v')W'j^by\B{v))  = ^{v  ■ v'){v  + vY 
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where  and  b^i  are  the  heavy  quark  fields.  The  heavy  quark  symmetry  used 
in  B physics  phenomenology  represents  substantial  progress  in  the  theoretical  de- 
scription of  B decays.  In  the  next  section  we  discuss  semileptonic  B decays,  for 
which  HQET  has  been  successfully  used  to  derive  decay  rates. 

1.3.3  Semileptonic  Decays  to  Mesons 

The  combination  of  large  branching  fractions  and  high  reconstruction  efficien- 
cies have  allowed  experiments  such  as  CLEO  to  measure  several  semileptonic  B 
decays  with  high  accuracy  [17].  This  wealth  of  experimental  results  has  allowed 
phenomenologists  to  compare  theory  and  experiment.  The  decay  kinematics  of  a 
specific  B semileptonic  decay  dictate  the  type  of  form  factor  contributions  to  the 
decay  rate.  In  the  case  of  .B  — D*W,  for  example,  there  are  no  corrections  to 
the  decay  rate.  This  behavior,  which  is  explained  by  Luke’s  theorem  [18],  implies 
that  the  HQET-derived  decay  rate  for  .B  — > D*li>  has  low  sensitivity  to  higher 
order  perturbative  corrections  as  well  as  non-preturbative  effects. 

1.3.4  Hadronic  Decays  to  Mesons 

Whereas  HQET  has  been  useful  in  describing  semileptonic  B decays,  an  equally 
accurate  description  of  hadronic  B decays  using  HQET  has  only  recently  begun  to 
be  pursued  for  two-body  decays  to  mesons  [19,  20,  21].  Whereas  in  semileptonic  B 
decays  one  of  the  two  currents  in  a matrix  element  is  weak,  and  therefore  calculable 
to  all  orders  in  perturbation  theory,  in  the  hadronic  case  we  have  matrix  elements 
of  four-quark  operators  with  hadronic  uncertainties  due  to  the  exchange  of  gluons 

and  quarks. 


20 


In  energetic  hadronic  two-body  B decays  hadronization  is  assumed  to  take  place 
after  the  quarks  that  form  each  of  the  two  hadrons  have  traveled  sufficiently  long 
distances  for  there  to  have  been  no  significant  exchange  of  gluons  or  quarks  between 
them.  This  decay  characterisitic  is  referred  to  as  factorization,  in  which  the  matrix 
elements  of  four-quark  operators  factorized  to  independent  current  elements  for 
each  hadron.  By  using  the  operator  product  expansion  (OPE)  [22,  23],  the  weak 
interaction  effects  are  treated  as  separate  from  the  long  range  strong  interaction 
effects.  The  HQET  effective  weak  hamiltonian  for  b c,u  transitions  after  this 

procedure  is  given  by: 

'^eff  = + C2(//)Q2^*’]  + h.C.)  + Qb-^u  + Qpenguin 

where  Ci,  2(/z)  are  scale  dependent  Wilson  coefficients.  Qb  — and  Qpenguin 

are,  respectively,  four-quark  operators  for  b ^ c,  b u,  and  penguin  decays  [20]. 
The  relative  strength  of  each  type  of  operator  as  well  as  the  Wilson  coefficients  are 
decay-dependent. 

Applying  the  factorization  hypothesis  to  He//  to,  for  example,  the  decay  am- 
plitude of  the  decay  D+tt",  results  in: 


A factorized  — (tJ'  ) 


where  Ci  can  be  verified  with  experiment.  As  the  number  of  hadronic  two-body 
B decays  and  the  accuracy  with  which  their  decay  rates  are  measured  increases 
in  the  near  future,  it  will  be  possible  to  test  the  decay  rates  derived  using  HQET 
similarly  to  how  it  has  been  done  for  semipletonic  B decays. 
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1.4  B Baryons 

A distinctive  feature  of  the  B meson  system  is  that  the  large  mass  of  the  b-quark 
allows  for  many  of  the  weak  decays  of  the  B meson  to  include  the  creation  of  a 
baryon-antibaryon  pair.  The  as  yet  unmeasured  decay  p n bars  this  feature 

from  being  unique.  In  Table  1.9  we  outline  the  B —¥  Baryons  decays  allowed  in 
the  Standard  Model. 


Table  1.9:  Representative  B Baryons  decay  mechanisms 


Quark-level  mechanism 

Sample  decay  mode 

Semileptonic  decays: 
b^cW-,W~ 
b-^uW~,W-  ^l~vi 
Hadronic  decays: 
b-^cW-,W-  ^cs 
b -4  cW~,  W~  -4  cs 
b-^cW~,W-  ud 

b-^cW-,W~  -4  ud 
b sg, 

b uW~ , W~  -4  ud 

B -4  Acpl~i'i 
B -4  ppl~vi 

B -4  J/ipK,  J/'ip  Baryons 
B :zcAc 
B AcPtt 
)■  D*~  p n 
B -4  Ap,B  — )■  ppK 
B — >■  pfin 

1.4.1  Results  to  Date 

The  B pX  was  measured  by  CLEO  to  be  8.0  ± 0.5  ± 0.3%[24],  assuming 
B -4  pX  = B nX.  Based  on  this  number  we  expect  roughly  8%  of  all  B 
mesons  in  our  data  to  be  a R -4  Baryons  event.  The  Ac,  Sc  and  :zc  charmed 
baryons  have  been  inclusively  measured  in  B decays  [25,  26].  Upper  limit  exclusive 
measurements  have  been  reported  to  date  for  B AcP^  with  Z = e, //  [27],  and 
selected  two-body  rare  decays  [28]. 
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Based  on  the  inclusive  measurements  reported  to  date,  B -4  Baryons  decays 
to  date  have  been  expected  to  be  produced  predominantly  in  decays  of  the  type 
B -4  AcpX,  and  this  is  the  type  of  decay  that  has  been  exclusively  reconstructed 
to  date  [29].  A typical  Feynman  diagram  for  a B ^ AcpX  decay  is  shown  in  Figure 

1.7. 


Figure  1.7:  A Feynman  diagram  for  jB  Kc^pt:  . 


1.4.2  The  Argument  for  B -4  [D]NNX  modes 

Combining  the  recently  measured  value  B{A'^  -4  pK  tt'*')  = 5.0  ± 0.5  ± 1.2% 
[30]  with  estimates  of  the  product  branching  fraction  B{B  Ac  A)  x B {Ac 
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pK~’n'^)  = 0.18%  [31]  we  can  determine  that  B AcpX  modes  account  for  only 
3.6%  of  5 — > Baryons,  which  is  approximately  half  of  the  total  B Baryons 
rate,  as  measured  in  the  inclusive  B -4  pX  measurement  [24]. 

Based  on  our  current  knowledge  of  J5  — ^ Baryons,  there  must  be  processes 
other  than  production  that  contribute  to  the  B —>■  Baryons  rate.  It  is  our  goal 
to  find  evidence  for  B — > Baryons  decays  which  do  not  involve  a Ag.  Dunietz 
[32]  suggested  that  modes  of  the  type  B -4  DNNX,  in  which  D represents  any 
charmed  meson,  and  N a proton  or  neutron,  are  likely  to  be  sizeable.  B —>  DNNX 
final  states  can  arise  either  from  the  hadronization  of  the  W boson  into  a baryon- 
antibaryon  pair,  or  from  the  production  of  a highly  excited  charmed  baryon  that 
decays  strongly  into  a baryon  plus  a charmed  meson.  CLEO  previously  reported 
an  inclusive  upper  limit  for  B{B  DNNX)  at  the  90%  confidence  level  of  < 4.8 
%[24]. 

1.4.3  Thesis  Overview 

In  this  thesis  we  will  attempt  the  exclusive  reconstruction  of  two  specific  B —> 
[D]NNX  decay  modes,  D*~  p p 7t+  and  D*~  p n.  Typical  Feynman 
diagrams  for  B^—^  D*~  p n and  4 D*~  p p tt"**  are  shown  in  Figures  1.8  and 

1.9,  respectively. 

The  choice  of  these  two  modes  is  guided  by  the  following  criteria: 

1.  D*  mesons  have  the  lowest  signal  contamination  among  the  D mesons. 

2.  Both  decays  can  occur  via  external  W decay.  Although  this  characteristic  is 
not  a principle,  to  date  only  the  b c decay  modes  that  have  been  measured 
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share  this  characteristic.  The  reasons  for  the  predominance  of  these  decays 
are  not  known. 

3.  These  are  the  two  modes  with  the  lowest  decay  daughter  multiplicity,  which 
translates  into  the  highest  reconstruction  efficiency. 

4.  jB°-4  D*~  p n and  D*~  p p have  low  combinatoric  backgrounds. 

We  report  here,  for  the  first  time,  evidence  for  decays  of  the  type  B -> 
[DjA^iVX,  and  present  measurements  of  the  branching  fractions  D*  p 

p TT"*")  and  B{B^-^  D*~  p fi).  The  charge  conjugate  process  is  implied  in  the  recon- 
struction of  B^-^  D*~  p p However,  in  the  reconstruction  of  > D*  p h, 
only  the  mode  with  the  antineutron  is  used  in  our  measurement  because  neutrons 
do  not  have  the  distinctive  annihilation  signature.  These  measurements  invalidate 
the  previous  assumption  that  B -4  Baryons  is  dominated  by  Ac  decays,  while  es- 
tablishing evidence  for  the  existence  of  a new  type  of  B decay  mechanism  with  a 

sizeable  decay  rate. 

The  thesis  is  divided  as  follows; 

In  Chapter  2 we  describe  the  CLEO  detector.  We  place  special  emphasis  on 
the  electormagnetic  calorimeter,  which  we  use  to  select  antineutron  candidates. 

In  Chapter  3 we  outline  the  selection  criteria  for  the  B^-^  D*~  p p 7t+  and 
D*~  p fi  decay  daughters.  We  place  special  emphasis  on  our  selection  of 

antineutron  showers  due  to  the  novelty  of  their  use. 

In  Chapter  4 we  describe  our  measurement  oi  B^-^  D*  p p tt"*"  , in  which  we 
use  a recontruction  technique  used  to  reconstruct  decays  in  which  the  energy  of  all 
the  decay  daughters  is  well  determined. 
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In  Chapter  5 we  describe  our  measurement  of  > D*~  p n , in  which  we 
use  a recontruction  technique  which  is  similar  to  that  used  to  reconstruct  other  B 
decays  with  missing  energy.  In  our  case  the  missing  energy  is  due  to  the  antineutron 
candidate. 

In  Chapter  6 we  conclude  by  summarizing  our  results,  stressing  on  their  sig- 
nificance, and  outlining  possible  decay  modes  that  we  believe  are  important  and 
measurable  with  the  expectedly  larger  datasets  available  to  future  studies. 


26 


Figure  1.8;  Two  Feynman  diagrams  for  D*  p n 
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Figure  1.9:  A Feynman  diagram  for  > D*  p p 


CHAPTER  2 
CLEO  II  DETECTOR 

The  CLEO  II  detector  [33]  is  centered  at  the  interaction  point  resulting  from 
the  collision  of  an  electron  beam  and  a positron  beam,  each  at  ~ 5.3  GeV  beam 
energy,  located  at  the  Cornell  Electron  Storage  Ring  (CESR).  The  beams  collide 
almost  head-on,  resulting  in  a center-of-mass  total  energy  of  ~ 10.6  GeV.  It  is  in 
the  energy  range  near  this  value  that  the  T resonances,  with  the  quark  content  bb, 
are  produced.  As  shown  in  Figure  2.1,  the  T(15)  to  T(4S')  resonances  are  found 
in  the  energy  range  9.44  GeV  to  10.62  GeV.  The  T(4«S')  is  above  the  threshold  for 
strong  decay  to  BB  pairs  to  take  place.  Close  to  100%  of  the  T(45)  decay  rate  is 

to  nearly  equal  numbers  of  charged  and  neutral  BB  pairs. 

Data  is  taken  at  the  T(45)  resonance  (to  be  referred  to  as  ON  resonance) 
to  study  B decays.  Because  approximately  2/3  of  the  ON  resonance  data  are 
composed  of  events  in  which  the  initial  qq  pair  is  not  an  T(4S'),  a sample  of  data 
is  also  taken  60  MeV  below  the  resonance  (to  be  referred  to  as  OFF  resonance)  m 
order  to  subtract  the  non-T(4S'),  or  continuum  component  of  the  ON  resonance 
data.  The  e+e'  annihilation  at  or  near  the  T(45)  resonance  yields  a wide  variety 

of  possible  final  states,  some  which  are  shown  in  Table  2.1. 

The  decay  product  of  an  e+e"  annihilation  is  called  an  event.  Because  final 

states  differ  greatly  in  cross  section  (the  frequency  with  which  events  of  a given 
topology  are  produced)  some  event  types  are  produced  with  high  frequency,  while 
other  event  types  are  produced  with  low  frequency.  Electronic  triggering  on  an 
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a (e^e'  Hadrons)(nb) 
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Figure  2.1;  Cross  section  into  hadrons  from  the  collision  of  e“*'e  beams  at  CESR 
as  measured  by  the  CLEO  II  detector  in  the  energy  range  9.44  GeV  to  10.62  GeV 
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Table  2.1;  Some  e+e“  annihilation  final  states 

e"''e“  -4  77 

e+e“  -)■  (with  l=e,  ^,r) 
e'^e~  ->  T(45)  — > BB  (requires  an  energy  threshold) 
e+e“  — > qq  (followed  by  hadronization) 


event-by-event  basis  is  used  to  select  some  or  all  of  a given  type  of  event.  Triggering 
allows  us  to  record  only  the  events  that  we  are  interested  in  studying,  most  of  which 
are  of  types  produced  with  low  cross  sections. 

2.1  Sub-detector  Components 

In  order  to  get  enough  information  out  of  an  event,  we  use  an  ensemble  of 
sub-detectors,  each  yielding  incomplete  information  about  the  event.  When  the 
information  from  all  sub-detectors  is  combined,  we  have  sufficient  information  to 
measure  useful  physics  properties.  A front  view  and  a side  view  of  the  CLEO  II 

detector  are  shown  in  Figure  2.2  and  Figure  2.3. 

From  the  innermost  to  the  outermost  (with  respect  to  the  beam  pipe,  which  is 

located  at  the  center  of  the  detector) , the  sub-detectors  are: 

1.  Vertex  detector 

(a)  PTL  (precision  tracking  layers)  detector,  used  during  the  earlier  part  of 
data  recording.  These  data  will  be  referred  to  as  CLEO  II  data. 

(b)  3-layer  SVX  (silicon  vertex  detector) , used  during  the  later  part  of  data 
recording.  These  data  will  be  referred  to  as  CLEO  II. 5 data. 


Figure  2.2:  Front  view  of  the  CLEO  II  detector 
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Magnetic  Coil 
Barrel  Crystals 
Time  of  Flight 

Endcap  Crystals 

Endcap  Time  of  Flight 


Figure  2.3i  Side  view  of  the  CLEO  II  detector 
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2.  Drift  chamber. 

3.  TOP  (Time-of-Flight)  detector. 

4.  Electromagnetic  calorimeter. 

5.  Muon  detector. 

The  volume  including  all  but  the  muon  detector  is  enclosed  in  a 1.5  Tesla 
superconducting  magnet.  An  important  feature  of  this  magnet  is  the  uniformity  of 
its  magnetic  field,  which  ensures  that  charged  particles  bend  uniformly  regardless 
of  where  in  the  detector  the  particle  travels.  A clear  introduction  to  detectors  as 
well  as  experimental  methods  in  high  energy  physics  is  found  in  Perkins  [34]. 

2.2  Tracking  System 

A charged  particle  traversing  a magnetic  field  in  the  presence  of  charged  wires 
in  a chamber  containing  gas  will  ionize  this  gas  as  it  loses  energy.  We  measure 
the  time  at  which  this  process  takes  place  as  well  as  the  energy  collected  by  each 
wire.  These  measurements  allow  us  to  know  the  position  of  the  particle  in  time 
and  the  energy  released  at  a number  of  points,  or  hits,  along  its  trajectory.  For  a 
given  momentum,  the  rate  at  which  a particle  loses  energy  along  this  trajectory, 
measured  as  dE/dx,  is  dependent  on  its  mass,  thus  allowing  us  to  separate  protons, 
kaons,  and  pions.  The  TOF  (Time-of-Flight)  of  a particle  in  a scintillating  medium 
is  also  dependent  on  its  mass  and  momentum.  TOF  measurements  yield  a second 
way  to  separate  protons,  kaons,  and  pions. 
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2.2.1  PTL  Detector 

The  PTL  detector  is  an  inner  drift  chamber  composed  of  six  layers  of  straw 
tubes.  There  are  64  axial  wires  for  each  layer,  and  there  is  a half  cell  stagger 
between  sequential  layers.  The  PTL  detector  does  not  measure  the  longitudinal, 
or  z-axis,  position  of  the  particle.  The  PTL  transverse  position  measurements  are 
more  precise  than  those  from  the  drift  chamber. 

2.2.2  SVX  Detector 

In  later  running  of  CLEO  II,  the  PTL  drift  chamber  was  replaced  by  a 3-layer 
SVX  detector  capable  of  longitudinal  as  well  as  axial  measurements  [35],  each 
measurement  taking  place  on  the  two  sides  of  each  of  the  silicon  wafers.  The 
radii  of  the  SVX  layers  are  2.35  cm,  3.25  cm,  and  4.75  cm,  for  layers  1,  2,  and  3, 
respectively.  It  is  composed  of  96  wafers  arranged  into  8 octants  of  12  wafers  each, 
with  26,208  data  readout  channels.  The  intrinsic  resolution  from  e+e“ 
events  at  normal  incidence  is  29  /^m. 

The  improved  measurement  resolution  of  the  SVX  detector  in  comparison  with 
the  PTL  detector  allows  for  more  accurate  determination  of  the  event  vertex.  This 
advantage  is  most  useful  to  lifetime  studies,  yet  it  does  not  affect  greatly  the  results 
presented  here. 

2.2.3  Drift  Chamber 

The  drift  chamber  system  (the  main  drift  chamber  and  the  vertex  detector), 
together  with  the  SVX  or  PTL,  are  used  to  measure  the  momentum  of  charged 
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particles.  Some  vertex  detector  and  drift  chamber  parameters  are  shown  in  Table 
2.2.  The  beam  pipe  is  located  at  radius  3.5  cm  in  the  CLEO  II  data  and  at  radius 
2.0  cm  in  the  CLEO  II. 5 data. 


Table  2.2:  Vertex  detector  and  drift  chamber  parameters 


Detector 

Layers 

Radius  (cm) 

Wires  per  layer 

PTL  (CLEO  II  only) 

6 

4.7  to  7.2 

64 

Vertex  Detector  (VD) 

10 

8.4  to  16 

64  (first  5),  96  (second  5) 

Outer  Drift  chamber 

51 

17.5  to  95 

96  to  384 

The  r — (j)  and  z measurement  resolutions  for  each  of  the  sections  are  shown  on 
Table  2.3. 


Table  2.3:  Drift  chamber  resolution 


Detector 

r — (f)  resolution 

z resolution 

PTL  (CLEO  II  only) 

90  jjm 

N/A 

Vertex  Detector  (VD) 

150  /im 

0.75  mm 

outer  drift  chamber 

110  ijm 

3 cm 

The  Vertex  Detector  (VD)  is  bounded  by  concentric  inner  and  outer  cathode 
strips  which  provide  z measurements.  The  segmentation  of  the  VD  cathode  strips 
is  5.85(6.85)  mm  along  z,  which  is  the  beam  direction,  on  the  inner(outer)  cath- 
ode. Segmented  cathodes  also  bound  layers  1 and  51  of  the  outer  drift  chamber. 
Segmentation  is  about  1 cm  along  z. 

f 

2.2.4  Momentum  and  Angular  Resolution 


There  are  two  factors  that  affect  the  track  momentum  resolution: 
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1.  The  error  in  the  measurement  of  the  track  curvature  due  to  the  hit-level 
measurement  error  in  drift  distance.  This  resolution  component  is  parametrized 
by  a term  linear  in  pt,  the  transverse  momentum. 

2.  Multiple  scattering  at  material  boundaries  which  cause  the  track  trajectory 
to  deviate  from  a helix.  This  resolution  component  is  parametrized  by  a constant 
term. 

The  parametrization  for  CLEO  II  data  is,  in  GeV: 

(^)2  ^ (0.0011pi)2  + (0.0067)2 

£^Z 

In  Table  2.4  we  show  this  resolution  in  MeV  for  selected  values  of  pt- 

Table  2.4:  Momentum  resolution  (5pt)  for  the  CLEO  II  data  at  selected  values  of 
Pt 


Pt 

ht 

0.5  GeV 

1.0  GeV 
1.5  GeV 

2.0  GeV 

3.3  MeV 
6.8  MeV 
10.4  MeV 
14.1  MeV 

The  angular  resolution  is  measured  using  e'^e"  — > pt^ pt~  events,  in  which  the 
typical  Pt  = 5.0  GeV.  8^—1  mrad  and  89  = A mrad.  We  expect  8(j)  and  86  to 
be  higher  for  the  tracks  we  use  in  our  analysis  since  muons  have  a much  lower 
probability  of  multiple  scattering  than  do  other  charged  particles. 

2.2.5  dE/dx  Measurements 

dE/dx  is  a function  of  particle  mass  and  momentum,  since  p — where 

P = v/c.  The  degree  of  separation  we  are  able  to  achieve  is  shown  in  Figure  2.4 
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for  the  CLEO  II  data.  The  51  layers  of  the  outer,  or  main,  drift  chamber— to  be 
referred  to  as  DR— are  used  to  measure  the  specific  ionization  energy  loss  (dE/dx) 
of  particles. 

A mixture  of  argon-ethane  gas  was  used  in  the  main  drift  chamber  for  the 
CLEO  II  data.  This  mixture  was  changed  to  helium-propane  during  data  taking 
for  CLEO  II. 5,  which  allowed  for  an  improvement  in  the  dE/dx  resolution,  resulting 
in  better  charged  particle  separation.  In  Figure  2.4  each  of  the  particle  bands  has 
been  plotted  after  a subtraction  of  some  of  the  higher  dE/dx  data  points  in  a given 
particle  band. 

The  main  drift  chamber  has  contiguous  cells  each  with  a sense  wire  surrounded 
by  field  wires,  as  shown  in  Figure  2.5.  Overall,  there  are  three  field  wires  for  every 
sense  wire  in  the  main  drift  chamber.  A number  of  corrections  are  applied  to 
optimize  the  resolution; 

1.  Dip  angle  saturation:  tracks  perpendicular  to  the  sense  wires  have  the  highest 
density  of  ionization  along  the  z direction.  The  amount  of  collected  charge 
is  reduced  by  electric  shielding  for  these  tracks. 

2.  Drift  distance:  varies  depending  on  the  field  configuration  of  each  cell. 

3.  (r,  (f>)  entrance  angle:  its  magnitude  as  well  as  its  sign  is  field  dependent. 

4.  Axial-stereo  layer:  cells  for  axial  and  stereo  layers  have  a different  field  de- 
pendence. 
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Momentum  (GeV/c) 


Figure  2.4: 


dE/dx  vs.  track  momentum 
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Figure  2.5:  Main  drift  chamber  (DR)  wire  arrangement 
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2.2.6  Time-of-Flight  Measurements 

Time-of-flight  detectors  surround  the  DR.  A bar  of  organic  scintillating  material 
5 cm  thick  has  photomultiplier  tubes  at  each  end  for  the  barrel  detectors,  and  at 
one  end  for  those  in  the  endcaps.  A time  measurement  is  made,  with  a 154  ps 
resolution  for  the  barrel,  and  272  ps  resolution  for  the  endcaps.  From  the  time  and 
distance  travelled,  a 1//?  quantity  is  defined.  1//3  varies  by  particle  type,  as  shown 
in  Figure  2.6. 

2.3  Electromagnetic  Calorimeter 

The  electromagnetic  calorimeter  is  used  to  measure  the  electromagnetic  energy 
deposition  of  charged  and  neutral  particles.  It  is  composed  of  7800  Csl  crystals.  A 
clustering  algorithm  is  used  to  combine  the  energy  deposition  in  a crystal  region, 
which  is  called  a shower. 

The  information  from  detector  components  is  used  in  our  analysis  in  a way 
that  is  consistent  with  previous  CLEO  II  analyses  with  the  exception  of  measure- 
ments from  the  electromagnetic  calorimeter.  The  electromagnetic  calorimeter  has 
been  used  previously  to  measure  electron  and  photon  energy  deposition.  In  our 
measurement  of  D*~  p n we  are  required  to  select  showers  that  are  consis- 
tent with  being  due  to  antineutrons  annihilating  with  the  Csl.  The  antineutron 
selection  procedure  is  successful  for  the  first  time  at  CLEO. 


TOF 
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Momentum  (GeV/c) 

Figure  2.6:  Time-of-Flight  vs.  track  momentum 
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2.3.1  Dimensions 

The  calorimeter  is  within  the  1.5  Tesla  magnetic  field.  All  crystal  faces  are 
at  1 m from  the  interaction  point,  facing  it  in  such  a way  that  showers  reach  all 
crystals  at  normal  incidence.  A partial  diagram  showing  some  of  the  barrel  and 
one  of  the  endcaps  is  shown  in  Figure  2.7.  Each  calorimeter  crystal  is  5-cm  x 

I 

5“Cm  X 30-cm,  where  the  latter  is  the  length  of  the  crystal.  The  choice  of  thalium 
doped  Csl  for  the  calorimeter  crystals  took  into  consideration  factors  such  as  cost, 
resistance  to  cracking,  high  density,  and  short  radiation  length  (1.83  cm).  Because 
the  calorimeter  is  16.4  radiation  lengths  deep,  1%  of  the  energy  of  a 5 GeV 
electron  leaks  out  of  it.  There  are  6,144  barrel  crystals  and  828  crystals  for  each 
endcap. 


MU 


1.74m 


1.44m 


1.02m 

0.95m 


BCC  Cables 
(_  BTF 


PTL,  VR, 
DR  Cables 

ECC  Cables 
ETF  Cables 


Figure  2.7:  Layout  of  CLEO  II  detector  showing  barrel  and  endcap  calorimeter 
sections 
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The  crystals  are  arranged  in  a rectilinear  grid  with  care  taken  to  have  a min- 
imum of  material  between  crystals.  The  material  in  front  of  the  barrel  section  is 
described  by  radiation  lengths  in  Table  2.5. 

Table  2.5:  Radiation  lengths  of  material  before  barrel  section  of  the  calorimeter 


Material 

Radiation  length  (%) 

beam  line  to  layer  1 of  DR 

2.5 

argon-ethane  gas  and  rest  of  DR 

0.4 

DR  outer  cathode  layer 

1.0 

outer  DR  wall 

2.0 

TOP  counters 

12.0 

There  is  more  material  in  front  of  the  endcap  calorimeter  sections,  which  de- 
grades measurement  quality.  We  use  endcap  information  in  the  7t°  reconstruction, 
but  not  for  the  antineutron  selection  criteria.  The  central  barrel,  which  covers  71% 
of  the  polar  angle  (45°  < ^ < 135°),  has  less  material  in  front  of  it,  thus  providing 
the  measurements  with  highest  resolution. 

2.3.2  Clustering 

The  algorithm  involved  in  associating  a group  of  nearby  cells  with  energy  above 
a threshold  in  the  calorimeter  is  called  clustering.  An  average  of  430  crystals  have 
energy  recorded  in  a hadronic  event  (either  BB  or  continuum).  The  raw  ADC 
count  measurements  need  to  be  calibrated  using  a crystal-to-crystal  calibration. 
The  sample  used  is  Bhabba  events  (e‘*‘e“  — > e'*'e“).  This  sample  has  high  statistics 
and  a beam  energy  constraint.  The  electrons  and  positrons  in  Bhabba  events 
deposit  almost  100%  of  their  energy  in  the  calorimeter.  Crystal  noise,  which  is 
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in  the  range  of  a few  MeV,  does  not  affect  this  sample  appreciably.  The  set  of 
constants  obtained  from  the  crystal-to-crystal  calibration  needs  to  be  updated 
only  every  few  months.  The  main  change  is  the  preponderance  of  a few  crystals  to 
become  noisy. 

A clustering  algorithm  for  CLEO  should  accomplish  the  following: 

1.  Match  tracks  to  showers. 

2.  Allow  7T°  reconstruction. 

3.  Allow  separation  of  photons  and  KlS  from  antineutrons. 

There  are  two  clustering  packages  in  the  CLEO  II  data:  CCFC,  and  XBAL. 
Both  accomplish  these  three  requirements.  Even  though  we  use  CCFC,  using 
XBAL  would  have  yielded  results  similar  at  the  few  % level.  Our  discussion  of 
clustering  is  limited  to  CCFC.  A discussion  of  clustering  using  XBAL  would  be 
similar  to  ours. 

In  order  to  separate  shower  clusters  from  noise,  only  cells  with  energy  > 10 
MeV  are  considered  as  candidates  for  the  center  of  a shower.  Only  cells  that  are 
at  most  two  cells  away  from  a cell  in  the  cluster  are  added  to  the  cluster,  which 
allows  cells  without  an  energy  measurement  to  be  located  inside  a cluster.  The 
number  of  cells  N that  is  considered  a cluster  is  a function  of  the  total  energy  of 
the  cluster.  For  example,  4 cells  correspond  to  a 25  MeV  shower,  17  cells  to  a 4 
GeV  shower. 

The  position  vector,  which  contains  the  directional  cosines,  is  determined  using 
a weighted  sum  of  the  energy  measured  with  each  cell,  using  the  geometric  center 
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of  each  to  define  a centroid.  Two  corrections,  one  lateral,  and  another  longitudi- 
nal, are  applied  to  this  position  vector.  These  corrections  account  for  calorimeter 
segmentation  as  well  as  depth  into  the  Csl  crystal  of  the  shower  center. 

The  cluster  energy  E and  angular  resolutions,  using  Monte  Carlo,  are: 

Barrel:  ^[%]  = + 1.9  — l.OE 

Endcap:  ^[%]  = + 2-5 

Barrel:  a(f,[mrad]  = -f  1.9 
Endcap:  a(j,[mrad]  = 

Barrel:  a$[mrad\  = 0.8a,psin{d) 

Endcap:  (7o[Tnrad]  = + 5.6 

Using  the  example  of  a photon  in  the  barrel,  we  list  in  Table  2.6  energy  and 
azimuthal  angular  resolutions  at  two  cluster  energies.  Their  low  values  are  a testa- 

4 

ment  of  our  ability  to  properly  reconstruct  electromagnetic  showers.  A by-product 
of  this  ability  is  our  success  in  separating  baryon-antibaryon  annihilation  showers 
from  electromagnetic  showers. 

Table  2.6:  Energy  and  angle  resolutions  for  a photon  in  the  barrel  at  two  values 
of  cluster  energy 


Cluster  energy 

f[%] 

a^[mrad] 

100  MeV 

3.8 

11 

5 GeV 

1.5 

3 
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2.4  Muon  Detector 

By  the  time  a particle  has  traversed  all  the  sub-detector  components  enclosed 
within  the  muon  detector,  all  particles  except  muons  have  deposited  most  of  their 
energy  in  the  calorimeter  and/or  have  been  deflected  by  the  magnetic  field  within 
the  drift  chamber.  Since  the  muon  detector  is,  above  everything  else,  a large 
volume  of  iron,  only  muons  are  expected  to  pass  through  a significant  section  of 
the  muon  detector. 

The  nuclear  absorption  length  \ for  a muon  in  iron  is  16.8  cm.  It  is  considerably 
less  for  all  other  charged  tracks  we  study.  The  maximum  depth  of  iron  to  be 
traversed  by  a particle  is  between  7.2Aj  and  lOAj.  We  do  not  use  the  muon  detector 
in  our  analysis. 


CHAPTER  3 
PARTICLE  SELECTION 


In  order  to  reconstruct  D*~  p p and  D*  p n we  need  to 

have  low  background  samples  of  all  the  decay  daughters.  Protons,  antiprotons  and 
pions  are  selected  as  single  tracks,  the  mD^-niDO  distribution  is  used  to  select  D*’s, 
and  shower  selection  criteria  is  used  to  select  antineutrons.  We  select  the  B decay 
daughters  from  ON  resonance  hadronic  events.  The  ON  resonance  candidates 
can  be  fully  or  partially  reconstructed  B mesons,  or  background  candidates  from 
ON  resonance  continuum  events.  The  OFF  resonance  candidates  can  only  be 
continuum  events  as  the  total  event  energy  is  below  the  energy  threshold  for  BB 
production  . We  use  the  MnJit  histogram  package  [36],  and  its  function  fitting 
utility,  MINUIT  [37],  both  of  which  are  widely  used  by  experimental  high  energy 
physicists. 

To  suppress  the  already  small  number  of  continuum  background  candidates  in 
our  reconstructions,  we  select  events  using  the  parameter  R2GL  [38].  R2GL  is 
a measure  of  how  the  momentum  is  distributed  for  the  event.  A high  value  of 
R2GL  corresponds  to  continuum  events,  in  which  the  initial  quarks  hadronize  to 
form  a two-body  decay.  The  BB  momentum  distributions  are  the  result  of  two 
separate  two-body  or  higher  decays  that  can  be  most  approximately  described  by 
a spherically  symmetric  momentum  distribution,  which  tends  to  yield  a low  value 
of  R2GL. 
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3.1  Data  Sample 

The  CLEO  data  is  collected  at  the  T(45)  resonance  (ON  resonance)  and  60 
MeV  below  the  resonance  (OFF  resonance).  Roughly  2/3  of  the  ON  resonance 
data  is  continuum  and  the  remaining  1/3  is  BB  pairs.  We  show  in  Table  3.1  the 
integrated  luminosity  by  dataset. 


Table  3.1:  The  CLEO  data 


Dataset 

ON  luminosity 

OFF  luminosity 

jBR’s 

CLEO  II 
CLEO  II.5 

3.14  fb~^  ' 
6.03  fb~^ 

1.61  fb~' 
2.94  fb~^ 

3.30  X 10^ 
6.35  X 10® 

3.2  Monte  Carlo  Sample 

We  use  Monte  Carlo  generated  events  to  model  the  data  in  order  to  study  the 
efficiency  of  our  selection  criteria  as  well  as  the  effect  of  backgrounds.  Two  types 

of  Monte  Carlo  samples  are  used: 

\ 

1.  Generic  Monte  Carlo,  in  which  many  decays  are  generated  using  our  current 
knowledge  of  B decays. 

2.  Signal  Monte  Carlo,  in  which  one  of  the  B mesons  in  the  event  decays  via  a 
predetermined  mode,  such  as  B^—^  D*  p n . The  remaining  B meson  in  the 
event  decays  according  to  the  prescription  used  to  generate  generic  Monte 

Carlo  events. 

The  generic  Monte  Carlo  sample  consists  of  B decays  of  type  b -)■  cW~.  A 
detailed  explanation  of  the  criteria  used  for  the  generic  Monte  Carlo  sample  is 
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found  in  the  Appendix  of  [24].  This  sample  does  not  include  any  B -)■  [D]NNX 
decays  because  until  now  we  have  had  no  evidence  for  their  existence.  The  B -4- 
Baryons  decays  in  the  generic  Monte  Carlo  sample  are  of  the  type  B — f U.W'^,  in 
which  n -4  Ac^iV,  where  N is  an  antiproton  or  an  antineutron.  The  conjugate 

modes  are  also  generated. 


3.3  Track  Selection 

A combination  of  measurements  is  used  to  select  quality  tracks  that  are  most 
likely  to  be  the  particle  type  we  need  for  a given  reconstruction-pion,  kaon,  or 
proton.  The  tracks  we  use  in  our  reconstruction  are  particles  that  are  stable  (in 
the  case  of  protons) , or  decay  outside  our  detector  (in  the  case  of  kaons  and  pions) . 
Therefore,  a track  is  assumed  to  begin  near  the  interaction  point,  which  is  within 
a solid  volume  defined  by  the  uncertainty  with  which  we  can  measure  the  beam 
spot.  This  volume  is  a few  microns  wide  in  x and  y and  a few  hundred  microns 

wide  in  z. 

3.3.1  Fitting  Algorithm 

Tracks  are  fitted  using  DR  hits  by  a Kalman  algorithm  that  is  based  on  the 
assumption  that  the  track  is  a helix  inside  a vacuum.  The  change  of  trajectory 
due  to  material  interaction  is  taken  into  consideration  in  the  track  fit.  The  fit  is 
performed  from  both  ends  of  the  track.  The  inward  fit  begins  a hit-to-hit  swim, 
adding  hits  to  form  a track,  analogous  to  the  way  beads  are  strung  to  form  a 
necklace.  Where  there  is  more  than  one  hit  to  a layer  of  the  drift  chamber,  the 
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hit  with  the  better  resolution  is  chosen.  The  outward  fit  begins  at  the  outermost 
layer  of  the  drift  chamber,  performing  a hit  swim  into  the  drift  chamber. 

3.3.2  Drift  Chamber  Track  Variables 

Variables  we  use  which  are  defined  with  main  drift  chamber  measurements 
(DR)  are: 

1.  PQCD:  signed  3-momentum  of  track.  This  vector  does  not  have  a correction 
assuming  the  mass  of  a particle  type. 

2.  ZOCD:  distance  of  closest  approach  along  the  z-axis.  The  z measurement  res- 
olution is  considerably  worse  than  the  r - (f>  measurement  resolution.  How- 
ever, ZOCD  is  effective  in  separating  tracks  from  the  interaction  point  from 
material  interaction  tracks  and  long  decay  particle  daughter  tracks. 

3.  DBCD:  radial  distance  of  point  of  closest  approach  to  the  interaction  point. 
DBCD  is  usually  more  accurately  determined  than  ZOCD. 

All  tracks  are  required  to  pass  the  following  track  quality  cuts: 

1.  PQCD  < 1 GeV:  DBCD  (in  meters)  < 0.005-(0.0038)  x \PQCD\. 

2.  Low  momentum  tracks  are  more  likely  to  have  a poorly  measured  DBCD. 
PQCD  > 1 GeV:  DBCD  (in  meters)  < 0.001. 


3.  ZOCD  < 0.05  meter. 


51 


3.3.3  The  TRKMNG  Package 

A software  package  is  used  to  reject  duplicated  tracks  from  curlers  [39].  A curler 
is  a low  momentum  track  which  is  bent  more  than  one  revolution  inside  the  drift 
chamber.  A curler  spirals  inside  the  drift  chamber  as  the  momentum  of  the  particle 
decreases.  The  tracking  algorithm  will  often  assign  a track  number  in  the  list  of 
tracks  for  each  half-revolution.  The  TRKMNG  package  selects  the  half-revolution 
that  is  most  likely  the  closest  to  where  the  track  began  to  curl,  rejecting  all  other 
half-revolutions.  The  variable  used  in  the  TRKMNG  Package  is  tng(track)  > 0. 
All  negative  numbers  of  tng(track)  are  for  the  redundant  half-revolutions. 

3.4  Particle  Separation 

By  population  size  alone,  all  tracks  are  pions.  The  purpose  of  particle  sepa- 
ration, or  particle  ID,  is  to  separate  tracks  that  are  likely  to  be  something  other 
than  a pion  from  the  overwhelmingly  more  numerous  pions.  Protons  are  produced 
with  a significantly  reduced  rate,  therefore  constituting  a small  background  to 
all  other  charged  particles.  The  predominant  backgrounds  to  protons  are  kaons, 
pions,  electrons  and  muons.  For  the  case  of  kaons,  pions  represent  the  largest  back- 
ground. The  kaon  in  has  more  phase  space  available  than  the  kaon  in 

and  D°^K‘^7r"7r+7r“,  which  allows  its  kinematic  separation  from 
backgrounds.  Even  for  the  low  Xp  in  D*’s  from  B decays,  the  kaon  in 
can  be  easily  separated  from  the  pion  by  the  decay  kinematics. 

Charged  particle  identification  is  accomplished  by  combining  the  specific  ion- 
ization (dE/dx)  measurements  from  the  drift  chamber  with  time-of-flight  (TOF) 
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measurements.  A normalized  probability  ratio  L{  is  used  for  charged  particle  sep- 
aration. 

Li  is  definded  as  Lj  = Pil{Ppion  + Pkam  + Pproton)',  Pi  is  the  particle  hy- 
pothesis Xi  probability  combining  dE/dx  and  TOF  measurements,  with  Xi  defined 

as: 

2 : {dE f dx)measured  — {dE f dx)exyected'\2  i f (TOF) measured  — {TOP) expected  12 

""  I C^dE/dx  •*  ^TOF  ■' 

o dE/dx, TOF  are  the  deviations  from  the  mean  and  vary  by  particle  type.  We 
define  the  respective  x?  for  each  particle  type  using  the  (TdE/dx,TOF  for  the  identical 
particle  type.  Lkaon,  for  example,  is  defined  using  the  adE/dx,TOFS  for  kaons.  In 
Table  3.2  we  show  how  our  criteria  varies,  depending  on  how  much  information 
is  available.  By  defaulting  to  a pion  we  are  setting  Pi  = 1.0.  Useful  dE/dx 
information  is  available  for  the  vast  majority  of  tracks  passing  the  track  quality 
requirements  outlined  in  section  3.3.2.  Low  momentum  particles  (<  200  MeV)  curl 
before  reaching  the  TOF  detectors. 

Table  3.2:  Particle  identification  criteria 


^ dE/dx 

OtOF 

Action 

< 3 

< 3 

use  both 

< 3 

> 3 

use  cr dE/dx 

> 3 

< 3 

use  otof 

> 3 

> 3 

default  to  pion 

We  place  loose  requirements  of  Lpion  > 0.001  on  all  pion  candidates,  and  the 

kaon  candidates  in  A'*' 7T~.  Kaon  candidates  in  A"*" tt  tt  , and  D — >A  tt  tt  tt 

are  required  to  satisfy  Lkaon  > 0.4.  Proton  candidates  in  A®— ^ D p p and 
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^ D*  p n are  required  to  satisfy  Lproton  > 0.9.  In  addition,  electron  back- 
grounds to  protons  and  antiprotons  are  rejected  using  the  variable  cut  R2ELEC 
< 0.  R2ELEC  is  a logarithmic  probability  that  a particle  is  an  electron.  Negative 
values  of  R2ELEC  are  least  likely  to  be  for  an  electron.  R2ELEC  is  defined  using  | 
from  the  electromagnetic  calorimeter,  which  is  very  close  to  1 for  electrons,  which 
leave  all  of  their  energy  in  the  calorimeter,  as  well  as  a non-negligible  number  of 

protons. 

3.5  TfO  Reconstruction 

The  7T°’s  used  to  reconstruct  D®— >R"^7r“7r°  decay  via  77.  The  tt^’s  in  this 
decay  are  soft,  with  momenta  < 500  MeV.  The  decay  angle  for  77  is  large 
enough  at  this  momentum  range  for  the  individual  7’s  not  to  overlap.  On  average 
there  are  7 tt^^^’s  in  a BB  event.  The  average  number  of  7t°’s,  as  of  yet  unmeasured, 
should  be  near  this  number.  Since  shower  multiplicity  increases  with  decreasing 
shower  energy,  as  tt®  momentum  decreases,  the  combinatoric  shower  background 
increases.  A minimum  momentum  of  100  MeV  is  required  of  our  7r°  candidates  to 
suppress  this  background. 

A shower  is  considered  a 7 if  the  ratio  between  the  summed  energy  of  a 3 by 
3 grid  of  shower  cells  over  the  same  are  extended  to  a 5 by  5 grid  of  shower  cells 
is  almost  1.  This  ratio  is  referred  to  as  E90E25.  Most  photons  have  most  of  their 
energy  concentrated  within  a 3 by  3 grid  of  shower  cells. 

We  require  at  least  one  of  the  7'^  fo  be  in  the  good  barrel  section  of  the 
calorimeter,  which  is  defined  as  |cos(0.y)|  < 0.71.  When  one  of  the  7 s is  not  in 
the  barrel,  we  require  that  it  not  be  from  the  overlap  region  between  the  barrel 


54 


and  endcap  sections  of  the  calorimeter  or  near  the  beam  pipe,  both  of  which  yield 
degraded  shower  energy  measurements.  We  also  require  the  candidate  7r°  mass  to 
be  within  2.5  o (tt®  ct  ~ 5 MeV)  from  135  MeV  central  value. 

3.6  D*  Reconstruction 

Finding  D*  candidates  with  high  efficiency  and  low  backgrounds  is  an  important 
part  of  this  analysis.  In  the  decay  D*  — )■  D^TTsoft  is  only  slightly  larger 

than  the  mass  of  the  tTj.  Since  the  resolution  is  smaller  than  the  mp» 

resolution,  we  use  the  former  to  select  D*  candidates.  We  verify  our  D*  yield  with 
a previous  analysis. 

In  order  to  reconstruct  B — > D*X  with  D*  -4  D^iTsoft  correctly,  we  check  our 
yield  for  the  scaled  momentum  Xp  range  0.25  < Xp  < 0.35.  It  is  in  this  range  that 
we  find  w 50  % of  D*’s  to  be  generated  in  signal  R°-4  D*~  p p 7t+  and  D*~ 
p n Monte  Carlo.  We  check  our  yield  with  a previous  CLEO  result  [40]. 

3.6.1  The  KNLIB  Fitting  Package 

The  KNLIB  fitting  package  is  a collection  of  kinematic  routines  used  for  track 
vertexing  [41].  We  use  it  to  vertex  Z)^’s  and  Z)*’s.  s are  vertexed  using  the 
KNLIB  routine  KNVTX3.  A new  invariant  mass  is  calculated  with  this  vertex,  and 
only  the  mass  range  1.70  to  2.03  GeV  is  allowed  to  pass  on  to  D*  reconstruction. 
To  form  a D*  vertex,  we  use  KNLIB  routine  KNBVXK  which  combines  the  D 
decay  products  and  the  tTj  track  to  make  a vertex-constrained  mass  from  which 
mp*-mpo  is  calculated.  We  use  KNBVXK  with  the  flag  IVOPT  set  to  1 for  the  tTj- 
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i.e.  the  tt*  is  constrained  to  pass  through  the  vertex  formed  by  the  other  particles, 
but  not  used  to  determine  that  vertex.  The  latter  method  offers  a noticeable 
improvement  over  other  options  tested,  including  use  of  track  when  vertexing, 
as  well  as  using  the  vertex  instead  of  the  individual  tracks  when  vertexing 
the  D*.  We  place  a loose  cut  of  < 100  vertex  fit  for  all  three  modes. 

3.6.2  Fit  Optimization 

We  reconstruct  D®’s,  selecting  candidates  near  the  mean  mass.  After  the  ad- 
dition of  a TTg,  we  select  candiates  within  a range  of  rri£)*-moo  values.  We  use  the 
1998  Review  of  Particle  Properties  values  for  mj)0,  1.8646  GeV,  and  mD*-m.£»o, 
0.1454  GeV,  as  our  central  values  [5].  The  double  Gaussian  fits  model 
better  than  the  single  Gaussian  fits. 

The  optimization  procedure  is  as  follows: 

1.  The  background  is  fitted  with  a second  order  Chebyshev  polynomial. 

2.  The  signal  is  fitted  with  a double  Gaussian  distribution. 

3.  The  continuum  subtracted  data  m^o  distribution  is  fitted  using  abs(m£>»- 
m£)O-0.1454)  < 0.002  GeV,  which  is  a wide  cut. 

4.  The  widths,  as  well  as  the  means  of  each  of  the  double  Gaussians  are  allowed 
to  float.  The  best  fits  are  obtained  with  non-zero  Ameajv’s.  ^mean  is  the 

* a 

difference  between  the  mean  value  of  each  of  the  two  Gaussians.  We  integrate 
the  double  Gaussian,  keeping  95%  of  the  signal  symmetrically  from  the  PDG 


mean. 
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5.  The  final  set  of  double  Gaussian  m£)0  and  rriD*-'n^D°  cuts  from  data  are  used 
in  the  B reconstructions,  and  the  same  cuts  are  applied  to  the  Monte  Carlo 
to  find  €mc,  the  signal  Monte  Carlo  reconstruction  efficiency. 

Shown  in  Figures  3.1,  3.2,  3.3,  3.4,  3.5  and  3.6  are  the  double  Gaussian  fits  for 
each  of  the  three  modes  for  either  dataset,  all  in  the  Xp  range  0.25  < Xp  < 

0.35. 
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MINUITx' Fit  to  Plot  210&3 

Att.  mdif  axis 

File;  */disl<1  /air/bnndx/dstar_cleo2.his 
Plot  Area  Total/Fit  5825.9  / 5825.9 
Func  Area  Total/Fit  5805.2  / 5805.2 

/=  11  4.8  for  140-  9 d.o.f., 

Errors  Parabolic 


25-JUL-99  16:38 
Fit  Status  3 
E.D.M.  7.769E-07 

C.L.=  84.1% 

Minos 


F unction 
NORM 
CHEB01 
CHEB02 
Function 
AREA 
MEAN 
SIGMA1 
AR2/AREA 
DELM 
SIG2/SIG1 


1 : Chebyshev  Polynomial  of  Order  2 
2.13491E+05  ±1.0441E+04 

0.20214  ±7.0776E-02 

-0.49237  ±9.5628E-02 

2:  Two  Goussians  (sigmo) 

4065.5  ±115.4 

0.14542  ±2.9342E-05 

1.36560E-04  ±3.4423E-05 


0.88944 

8.63015E-05 

3.5746 


±3.9038E-02 

±3.4756E-05 

±0.7962 


O.OOOOE+00 

O.OOOOE-fOO 

O.OOOOE+00 

O.OOOOE+00 

Q.OOOOE+00 

O.OOOOE-t-00 

■O.OOOOE+00 

O.OOOOE+00 

■0.0000E±00 


+ O.OOOOE+00 
+ 0.0000E+00 
+ 0.0000E+00 

+ 0.0000E+00 
+ 0.0000E+00 
+ 0.0000E+00 
+ 0.0000E+00 
+ 0.0000E+00 
±0.0000E+00 


Figure  3.1*.  in  GeV  for  B — Z)*A  with  in  CLEO  II 
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MINUITx' Fit  to  Plot  210&3 

Atv.  mdif  axis 

File*.  */disk1  /air/bnndx/dstar_c!eo25.his 
Plot  Area  Totol/Fit  6783.2  / 6783.2 
Func  Area  Totol/Fit  6770.0  / 6770.0 


/=  126 

Errors 

Function  1 ! 

NORM 

CHEB01 

CHEB02 

Function  2! 

AREA 

MEAN 

S1GMA1 

AR2/AREA 

DELM 

S1G2/S1G1 


.8  for  1 40  — 9 d.o.f., 

Parabolic 

Chebyshev  Polynomial  of  Order  2 
2.10956E+05  ±1.2915E+04 

0.29909  ±6.6768E-02 

-0.42591  ±0.1134 

Two  Goussions  (sigmo) 

5083.6  ±142.3 

0.14553  ±1.0554E-05 

2.26699E-04  ±1.6435E-05 


0.55532 

3.45832E-05 

2.5902 


±5.21 12E-02 
±2.8669E-05 
±0.1501 


26-JUL-99  10:09 
Fit  Status  3 
E.D.M.  2.662E-07 

C.L.=  58.7% 

Minos 

-O.OOOOE+00  +0.0000E+00 
-O.OOOOE+00  +0.0000E+00 
-O.OOOOE-tOO  +0.0000E+00 

-O.OOOOE+00  +0.0000E+00 
-O.OOOOE+00  +0-0000E+00 
-O.OOOOE+00  +O.OOOOE+00 
-O.OOOOE+00  +0.0000E+00 
-O.OOOOE+00  +0.0O0OE+00 
-O.OOOOE+00  +0.0000E+00 


Figure  3.2:  tyid^-ttido  in  GeV  for  B D*X  with  in  CLEO  II. 5 
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MINUIT  Fit  to  Plot  210Sc3 

An.  mdif  axis 

File!  */disk1  /air/bnndx/dstar_cleo2.his 
Plot  Area  Total/Fit  24953.  / 24953. 

Func  Area  Total/Fit  24954.  / 24954. 

v^=  100.2  for  100-  9 d.o.f., 

Errors  Porabolic 

Function  1 ! Chebyshev  Polynomial  of  Order  2 
NORM  1.68800E+06  ±4.0965E+04 

CHEB01  0.42131  ±1.8368E-02 

CHEB02  -0.21635  ±3.6013E-02 

Function  2!  Two  Goussions  (sigma) 

AREA  6856.6  ±555.8 

MEAN  0.14552  ±3.5495E-05 

SIGMA1  3.69304E-04  ±6.8050E-05 

AR2/AREA  0.56666  ±0.1691 

DEEM  -6.31264E-05  ±9.0940E-05 

SIG2/SIG1  2.3621  ±0.3285 


27-JUL-99  10103 
Fit  Status  3 
E.D.M.  1.348E-06 

C.L.=  24.0% 

Minos 

-O.OOOOE+00  +0.0000E+00 
- O.OOOOE+00  + O.OOOOE+00 
-O.OOOOE+00  +O.0O00E+00 

-O.OOOOE+00  +0.0000E+00 
-O.OOOOE+00  +0.0000E+00 
-O.OOOOE+00  +0.0O00E+00 
-O.OOOOE+00  +0.0000E+00 
-O.OOOOE-hOO  +0.0000E+00 
-O.OOOOE+00  +0.0000E+00 


Figure  3.3;  mp.-mDO  in  GeV  for  B -)■  D*X  with  D°^A:+7r  7T°  in  CLEO  II 


j 
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MINUITx^  Fit  to  Plot  210Sc3 

An.  mdif  axis 

File!  */disk1  /air/bnndx/dstar_cleo25.his 
Plot  Area  Total/Fit  28444.  / 28444. 

Func  Area  Total/Fit  28432.  / 28432. 


Errors 
Function 
NORM 
CHEB01 
CHEB02 
Function 
AREA 
MEAN 
SIGMA1 
AR2/AREA 
DELM 
SIG2/SIG1 


70.7  for  100-9  d.o.f.. 

Parabolic 

1 ; Chebyshev  Polynomial  of  Order  2 
1.93601E+06  ±3.7300E+04 

0.42646  ±1.7573E-02 

-0.23873  ±2.7967E-02 

2\  Two  Goussions  (sigma) 

7531.3  ±454.4 

0.14552  ±1.7005E-05 

3.38045E-04  ±2.681  IE-05 


0.46416 

1.43254E-04 

2.7650 


±6.5637E-02 

±7.7580E-05 

±0.3587 


27-JUL-99  10147 


E.D.I 

C. 

Minos 

-O.OOOOE+00 
-O.OOOOE+00 
-O.OOOOE+00 

-O.OOOOE+00 
-O.OOOOE+00 
-O.OOOOE-fOO 
-O.OOOOE+00 
- O.OOOOE+OO 
-O.OOOOE+00 


Fit  Status 
h.  2.196E-06 

L.=  94.3% 


O.OOOOE+00 
O.OOOOE+00 
+ 0.0000E+00 

+ 0.0000E+00 
+ 0.0000E+00 
+ 0.0000E+00 
+ 0.0000E+00 
+ 0.0000E+00 
+ O.OOOOE+00 


Figure  3.4:  mo. -moo  in 


GeV  for  B -)■  D*X  with  D°->Ar+7r  7t°  in  CLEO  II.5 
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MINUITx' nt  to  Plot  210Sc3 

An.  mdif  axis 

File!  »/disk1  /air/bnndx/dstar_cleo2.his 
Plot  Area  Total/Fit  15419./  15419. 
Func  Areo  Total/Fit  15397./  15397. 


v'=  101.5  for  100-9  d.o.f., 

Errors  Parabolic 


Function  1 ! Chebyshev  Polynomiol  of  Order  2 
NORM  1. 06391 E+06  ±1.9036E+04 

CHEB01  0.38063  ±2.2132E-02 

CHEB02  -0.16830  ±2.9278E-02 


Function  2! 

AREA 

MEAN 

SIGMA1 

AR2/AREA 

DEEM 

SIG2/SIG1 


Two  Gaussions  (sigma) 


4161.1 
0.14550 

3.1  1927E-04 
0.51012 
5.25295E-05 
2.2124 


±196.5 

±2.0141E-05 

±4.2942E-05 

±0.1495 

±7.4998E-05 

±0.2508 


27-JUL-99  13105 
Fit  Status  3 
E.D.M.  1.064E-06 


C.L.=  21.2% 


Minos 


O.OOOOE-fOO 

O.OOOOE400 

O.OOOOE+00 

O.OOOOE-fOO 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE400 

O.OOOOE+00 

0.00O0E4O0 


+ 0.0000E+00 
+ 0.0000E+00 
+ O.OOOOE+00 

+ 0.0000E+00 
±0.0000E±00 
-fO.OOOOE+00 
+ 0.0000E+00 
+ 0.0000E±00 
4 0.0000E+00 


0.140  0.1425  0.145  0.1475  0.150 


Figure  3.5:  in  GeV  for  5 — > D *X  with  tt+tt  in  CLEO  II 
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MINUIT  x' Fit  to  Plot  210&3 

Att.  mdif  axis 

File!  */disk1  /oir/bnndx/dstar_cleo25.his 
Plot  Area  Total/Fit  1 0721 ./  1 0721 . 

Func  Area  Total/Fit  10704.  / 10704. 


27-JUL-99  13:47 
Fit  Status  3 
E.D.M.  3.768E-06 


F=  68.4  for  100- 

rrors 


9 d.o.f., 

Parabolic 


C.L.=  96.3% 

Minos 


Function  1 ! Chebyshev  Polynomial  of  Order  2 
NORM  6.86301  E+05  ± 1 .3867E-I-04  - O.OOOOE+00 

CHEB01  0.41991  ±2.7363E-02  - O.OOOOE+OO 

CHEB02  -0.22118  ±3.2620E-02  -O.OOOOE+OO 


Function  2!  Two  Gaussians  (sigma) 
AREA 
MEAN 
S1GMA1 
AR2/AREA 
DEEM 
SIG2/S1G1 


3335.0 
0.14552 
1.87938E-04 
0.70392 
9.37545E-06 
2.3674 


± 1 14.4 

±2.0129E-05 

±3.1388E-05 

±0.1019 

±3.8689E-05 

±0.2952 


- O.OOOOE+OO 
-O.OOOOE+OO 
-O.OOOOE+OO 
-O.OOOOE+OO 
-O.OOOOE+OO 
-O.OOOOE+OO 


+ O.OOOOE+OO 
+ O.OOOOE+OO 
+ O.OOOOE+OO 

+ O.OOOOE+OO 
+ O.OOOOE+OO 
+ O.OOOOE+OO 
+ O.OOOOE+OO 
+ O.OOOOE+OO 
+ O.OOOOE+OO 


Figure  3.6:  in  GeV  for  B — > D*X  with  tt^tt  in  CLEO  II.5 
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Shown  in  Table  3.3  are  the  rriDO  and  mD*-m£)0  double  Gaussian  data  cuts  for 
all  three  modes  and  for  each  dataset. 

Table  3.3:  rnoo  and  double  Gaussian  data  cuts 


Mode 

Dataset 

m£»o 

in  ± MeV 

TRD*~Tn£}0 

in  ± MeV 

CLEG  II 

17.5 

1.15 

CLEG  II.5 

15.0 

1.10 

CLEG  II 

26.0 

1.50 

CLEG  II.5 

27.0 

1.50 

7T'*'7r~ 

CLEG  II 

14.0 

1.30 

.D  ° A 7T  - 7T  7T  “ 

CLEG  II.5 

12.0 

0.90 

3.6.3  Comparison  with  B — > D*X 

To  verify  the  accuracy  of  our  D*  reconstruction,  we  compare  our  results  with 
the  PDG  value  Bt{B  D*X)  = 22.7  % [5].  Since  €mc  varies  with  momentum, 
we  need  to  compare  our  fully  reconstructed  sample  with  the  expected  number  in 
the  same  momentum  range.  In  B^^  D*~  p n and  B^^  D*~  p p signal  Monte 
Carlo  is  ~ 50  % of  the  generated  D*’s  are  in  the  scaled  momentum  (Xp)  range 
from  0.25  to  0.35.  Approximately  32  % of  Bx{B  D*X)  lies  in  this  momentum 
range  [40].  The  expected  D*  yield  is  given  by: 

(Number  of  charged/neutral  B’s)  x 

(Br(B  D*  X)  in  Xp  range  0.25  to  0.35)  x 

(Br(D°  mode))  x 

(Br(T»*  ^ D^Tisoft)  = 68.3  %)  x 

{D*  €mc)  X 
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(0.32/0.50:  efficiency  correction). 


This  number  is  efficiency  corrected  for  single  and  double  Gaussian  signal  fits 

using  D*~  p n Monte  Carlo.  The  corrected  product  is  the  expected  yield 

that  we  compare  to  our  respective  fitted  yield  for  each  mode. 

The  close  agreement  between  the  expected  yield  from  Monte  Carlo  and  the 

data  result  for  each  mode,  as  shown  on  Tables  3.4,  3.5,  and  3.6,  is  a function  of 

how  well  our  Monte  Carlo  models  D*  decays  from  B mesons  and  of  the  accuracy 

of  our  reconstruction  code.  The  double  Gaussian  fits  yield  better  agreement  than 

the  single  Gaussian  fits.  There  is  a significant  drop  in  reconstruction  efficiency  for 

all  three  modes  in  the  CLEO  II. 5 data.  This  drop  is  due  to  the  significantly 

reduced  reconstruction  efficiency  of  the  soft  pion  in  D*  — > D^iisoft  in  CLEO  II. 5. 

Whereas  in  CLEO  II  the  PTL  drift  chamber  is  used  in  the  tracking  algorithm, 

the  SVX  silicon  detector  which  replaces  it  in  the  CLEO  II. 5 data,  in  addition  to 

« 

having  more  material,  is  not  used  in  the  tracking  algorithm. 


Table  3.4:  D*  yield:  D^-^K^'k 


Dataset 

Single 

Gaussian 

^MC 

Expected 

yield 

Found 

Double 

Gaussian 

^MC 

Expected 

yield 

Found 

CLEO  II 
CLEO  II.5 

_ 30.8  % 
' 16.3  % 

3,862.3 

3,933.1 

4,117.8 

4,455.2 

33.0  % 
20.7  % 

4,138,1 

4,994,9 

4.065.5 

5.083.6 
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Table  3.5:  D*  yield:  7t° 


Dataset 

Single 

Gaussian 

Expected 

Found 

Double 

Gaussian 

Expected 

Found 

s/ 

^MC 

yield 

^MC 

yield 

VS  + B 

CLEO  II 

14.0  % 

6,371,3 

5,955.8 

15.3  % 

6,963.0 

6,856.6 

59.7 

CLEO  II.5 

8.4  % 

7,356.0 

6,138.5 

8.9  % 

7,793.9 

7,531.3 

61.4 

Table  3.6:  D*  yield:  T>0-^i^+7r-7r+7r- 


Dataset 

Single 

Gaussian 

^MC 

Expected 

yield 

Found 

Double 

Gaussian 

€-MC 

Expected 

yield 

Found 

s/ 

VsTb 

CLEO  II 
CLEO  II.5 

14.5  % 

6.5  % 

3,560.6 

3,071.3 

3,637.6 

3,100.5 

16.0  % 
6.9  % 

3,928.9 

3,260.3 

4,161.6 

3,335.0 

48.3 

47.7 

3.7  Antineutron  Showers 

We  need  to  define  a set  of  criteria  which  allows  us  to  select  antineutrons  with 
high  accuracy  without  incurring  a large  loss  in  efficiency.  The  following  character- 
istics, some  limiting,  some  exploitable,  apply  to  antineutron  showers  in  the  CLEO 
II  electromagnetic  calorimeter: 

1.  The  antineutron  shower  yields  an  incomplete  measurement  of  its  energy  and 
an  accurate  measurement  of  its  direction.  We  use  the  antineutron  shower 
energy  to  select  candidates,  but  cannot  use  this  energy  when  reconstructing 
the  B candidate.  We  can  use  the  well-measured  shower  direction. 
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2.  Antineutrons  frequently  annihilate  with  matter  in  the  calorimeter.  Since 
antineutron  annihilation  showers  have  distinctive  characteristics  that  enable 
us  to  separate  them  from  other  showers,  we  use  these  characteristics  to  select 
them. 

3.  By  baryon  number  conservation,  if  the  event  has  an  antineutron,  it  must  have 
a corresponding  baryon.  In  an  exclusive  reconstruction,  as  is  ours,  the  selec- 
tion of  an  proton  increases  the  probability  of  there  being  an  antineutron  in 
the  event  substantially  due  to  baryon  number  conservation.  In  addition,  once 
we  have  selected  a D*,  only  shower  candidates  within  a narrow  momentum 
cone  can  be  selected  in  the  event. 

We  refer  the  reader  to  other  studies  of  baryon-antibaryon  annihilation  [43,  44]. 

In  Table  3.7  we  outline  the  types  of  showers  encountered  in  a hadronic  event. 


Table  3.7:  Shower  types  and  energy  measured  in  calorimeter 


Particle 

Shower  type 

Energy  measured 

7T=^,  K^,  p+ 
e±,7 

Kl 

neutron 
p and  n 

minimum  ionizing 
electromagnetic 
soft  annihilation 
small  electromagnetic 
medium  to  hard  annihilation 

small  fraction 
full  measurement 
small  fraction 
very  small  fraction 
small  fraction 

We  are  limited  by  the  absence  of  an  independent  antineutron  sample  that  we 
can  study  to  define  our  selection  criteria.  However,  antipfotons  also  annihilate  with 
nucleons  in  the  calorimeter.  Therefore,  antiproton  annihilation  showers  compose 
the  shower  sample  in  data  and  Monte  Carlo  which  we  use  to  define  our  antineutron 
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selection  criteria  as  well  as  to  gauge  how  well  the  Monte  Carlo  models  antineutron 
annihilation  showers. 

3.7.1  Shower  Parameters 

We  use  the  CCFC  clustering  package  to  select  antineutron  candidates.  The 
CCFC  shower  package  was  not  optimized  to  separate  annihilation  showers  from 
other  showers  in  the  calorimeter.  Nevertheless,  we  find  that  the  parameters  previ- 
ously optimized  to  include  photons,  such  as  E90E25,  the  list  of  nearby  showers, 
and  track-to-shower  matching,  are  useful  in  excluding  them  as  well. 

The  shower  parameters  used  are: 

1.  E90E25  < cut  1,  a pre-determined  value:  E90E25  has  been  defined  in  our 
discussion  of  7t°’s.  It  is  very  close  to  1 for  photons,  and,  we  find,  farthest  from 
1 for  annihilation  showers,  overlapped  7r°’s-which  are  not  merged,  but  have 
the  two  photon  showers  very  near  each  other,  and  K^s.  An  energy  dependent 
cut,  cut  1,  on  E90E25,  is  applied  to  reject  99%  of  isolated  photons. 

2.  NNESH:  List  of  nearby  showers  (which  does  not  include  the  shower  to  which 
the  list  belongs).  The  area  encompassed  increases  with  energy.  It  is  in  this 
list  that  we  find  showers  near  the  antineutron  that  are  most  likely  to  be 
hadronic  split-offs.  We  call  the  shower  associated  with  the  list  the  main 
shower,  and  sum  the  energy  of  it  and  all  the  others  in  NNESH,  we  call  a 
group.  In  5°-)-  D*~  p n CLEO  II  signal  MC,  showers  in  the  NNESH  list 
are  tagged  to  the  parent  shower  in  excess  of  93  % of  the  time. 
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3.  |(cos(^)|:  The  calorimeter  can  be  divided  by  polar  angle  9 into  sections:  good 
barrel,  bad  barrel,  barrel/endcap  overlap,  bad  endcap,  good  endcap,  and  near 
beampipe.  We  use  the  good  barrel  section  throughout,  which  corresponds  to 

I(cos(0)l  < 0.71. 

4.  Match  level:  The  track  to  shower  match  levels  in  CCFC’s  NTRSH  array  are: 

(a)  Match  level  1:  shower  center  < 8 cm  from  the  track  projection. 

(b)  Match  level  2:  not  level  1,  but  > 1 member  crystal 

(c)  Match  level  3:  not  level  1 nor  level  2. 

We  use  match  3 showers  for  antineutron  candidates  and  match  1 and  2 show- 
ers for  antiprotons. 

3.7.2  Antiproton  Showers  in  Data 

Lacking  a sample  of  antineutrons  to  study  in  data  which  is  independent  of 
the  sample  we  will  use  to  measure  D*~  p n and  J5°— D*  p p tt'*'  , we  use 
antiproton  annihilation  showers  which  have  been  matched  to  a track.  Antiprotons, 
like  antineutrons,  annihilate  with  matter  in  the  calorimeter.  The  annihilation 
group  of  showers  that  results  from  this  process  has  a distinctive  signature:  the 
antibaryon  interaction  is  associated  with  a main  shower,  which  has  most  of  the 
energy  of  the  annihilation,  and  the  hadronic  splitoffs  from  the  interaction  are 
associated  with  satellite  showers  near  the  point  of  interaction,  each  making  a small 
contribution  to  the  group  energy. 

In  Figure  3.7  we  plot  shower  energy  Emam  vs.  signed  momentum  (PQCD)  for 
protons  and  antiprotons,  allowing  all  values  of  E90E25,  in  the  CLEO  II  data  from 
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a D*  skim.  Continuum  and  BB  hadronic  events  are  combined  in  this  plot.  The 
use  of  a D*  skim  is  a result  of  convenience. 

The  selection  criteria  used  to  generate  Figure  3.7  is  as  follows: 


1.  \PQC D{track)\  > 300  MeV. 


2.  Lproton  ^ 0.9. 


3.  Track-to-shower  match  level  1 


4.  dE / dx^eiectvon  ^ 2.2  to  reject  electron  fakes. 


5.  |(cos(0)|  < 0.71:  good  barrel. 


6.  All  values  of  E90E25  allowed. 


In  Figure  3.8  we  use  the  same  selection  criteria  as  was  used  for  Figure  3.7  with 
the  exception  of  a rejection  cut,  E90E25  < cut  1,  to  suppress  non-annihilation 
showers.  The  antiproton  annihilation  showers  are  typically  those  with  Emain  > 
300  MeV  and  in  the  track  momentum  range  500  MeV  to  900  MeV.  Note  that  for 
these  showers  there  is  only  a loose  correlation  between  the  reconstructed  energy  of 
the  shower  and  the  momentum  of  the  antiproton  that  produced  it.  We  can  identify 
an  annihilation  shower,  but  have  a very  poor  measurement  of  the  momentum  of 
the  particle  that  produced  it. 

The  horizontal  line  at  shower  energy  200  MeV,  sloping  upwards  at  PQCD  1 
GeV,  corresponds  to  minimum  ionizing  protons  and  antiprotons.  The  rest  of  this 
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File!  ♦/diskl /olr/dnnclx/onprot^. cieo2.his 

l[J  IDB  SyrntJ  Dole/Time  Areo  Mean  R.M.S. 

782  1 1 99C7i0/1023  2,835.^E+04  - 1 .0707E-02  1.026 

0.6473  0.4285 


Figure  3.7:  Emain  vs.  PQCD  for  protons  and  antiprotons  in  CLEO  II 

line  in  the  lower  momentum  range  corresponds  to  minimum  ionizing  pions.  In  the 
case  of  protons,  before  an  E90E25  cut  is  applied,  the  diagonal  line  corresponds  to 
protons  captured  in  the  calorimeter.  1.0  GeV  momentum  is  the  threshold  for  this 
capture.  Neutrons,  analogously  to  protons,  do  not  annihilate,  and,  since  we  have 
neither  momentum  information  nor  distinguishable  showers,  the  study  of  B decays 
with  neutrons  is  beyond  CLEO’s  capabilities. 


3.7.3  Antineutron  Selection  Criteria 

Using  antiproton  annihilation  showers,  we  devise  the  antineutron  selection  cri- 
teria to  be  used  when  reconstructing  D*~  p h , as  shown  in  Table  3.8. 

Although  this  selection  criteria  allows  a shower  to  have  no  nearby  daughters,  in 
which  case  Emain  = E group,  the  majority  of  antineutron  showers  have  at  least  one 
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rile;  ‘/disk l/oir/bnn(Jx/oopfot_uleo2. his 

ID  IDB  Symb  Dotc/Time  Areo  Meon 

782  1 1 990730/1023  1.3408E+04  -0.4700 

0.8003 


R.M.S. 

0.9499 

0.4182 


Figure  3.8;  Emain  vs.  PQCD  for  protons  and  antiprotons  in  CLEO  II 

nearby  daughter  in  Monte  Carlo,  as  well  as  in  our  exclusively  reconstructed 
signal  events. 

Table  3.8:  Antineutron  shower  selection  criteria. 


Track-to-shower  match  level  3 
E90E25  < cut  1 
jcOs(^s/i,o^er)  1 0.71 

Emain  > 500  MeV 
Egroup  > 800  MeV 


We  test  this  selection  criteria  in  a generic  BB  Monte  Carlo  before  and  after  the 
requirement  that  there  be  a proton  with  \PQCD{track)\  > 300  MeV  and  Lproton 
> 0.9  in  the  event.  The  generic  Monte  Carlo  sample  is  discussed  in  Section  3.2.  In 
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Table  3.9  we  show  the  results  before  and  after  applying  the  antineutron  selection 
criteria  cuts  for  a sample  before  the  proton  requirement  is  applied. 


Table  3.9:  Shower  population  in  a generic  BB  Monte  Carlo  sample  without  proton 
requirement 


Shower 

type 

Annihilation  cuts 
without  Egroup  cut 

Annihilation  cuts 
with  Egroup  > 800  MeV 

7 from  7T*^ 

44.5 

46.2 

Kl 

33.0 

22.2 

P 

0.5 

0.6 

1 

K- 

3.3 

1.9 

other(cc>,etc) 

1.4 

1.3 

ft 

17.3 

27.8 

TOTAL 

100.0 

100.0 

According  to  our  generic  Monte  Carlo  simulation,  even  before  a proton  require- 
ment is  applied  27.8%  of  the  annihilation-like  showers  in  a BB  event  are  antineu- 
trons. We  next  apply  the  proton  requirement  to  a generic  BB  Monte  Carlo  sample. 
The  results  are  shown  in  Table  3.10. 

Table  3.10:  Shower  population  in  a generic  BB  Monte  Carlo  sample  with  proton 
requirement 


Shower 

type 

Annihilation  cuts 

^Vlthout  group  cut 

Annihilation  cuts 
with  Egroup  > 800  MeV 

7 from 

16.3 

11.4 

Kl 

9.5 

3.9 

P 

2.7 

1.8 

2.0 

0.9 

other  (cc, etc) 

0.7 

0.5 

n 

68.8 

81.5 

TOTAL 

100.0 

100.0 
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The  above  results  are  not  surprising:  use  of  baryon  number  conservation  by 
selecting  a proton  increases  the  probability  of  an  event  having  an  antineutron  from 
a few  % to  50  %,  while  leaving  all  other  annihilation  like  backgrounds  at  near  the 
same  level.  Use  of  Egroup  cut  reduces  Kl  contribution  by  1/3  in  the  Non-baryon 
sample  and  by- 1 /2  in  the  baryon  sample. 

3.7.4  Antineutron  Backgrounds 

In  Figures  3.9  and  3.10  we  compare  the  shower  energy  spectrum  of  each  of  the 
two  major  backgrounds  to  antineutrons-7  from  tt'^’s,  and  Ki~to  that  of  antineu- 
trons. All  antineutron  selection  criteria  cuts,  including  Egroup  > 800  MeV,  have 
been  applied.  Figure  3.9  is  for  case  before  proton  requirement  has  been  applied, 
and  Figure  3.10  is  for  the  case  after  proton  requirement  has  been  applied.  The 
sample  is,  again,  generic  BB  Monte  Carlo.  The  backgrounds  that  have  not  been 
included  make  up  (4.3,3.0)  % of  the  (before,after)  proton  requirement  distribu- 
tions. In  Figure  3.9  the  solid  distribution  is  for  n’s,  the  slanted  lines  distribution 
is  for  7’s  from  tt^^’s,  and  the  white  distribution  is  for  KEs. 

In  Figure  3.10  the  white  distribution  is  for  h’s,  the  solid  distribution  is  for  7’s 
from  7T°’s,  and  the  slanted  lines  distribution  is  for  KEs.  The  large  contamination 
in  the  low  range  of  the  Emain  spectrum,  as  shown  in  Figure  3.9  and  Figure  3.10 
does  not  affect  D*~  p n since  for  this  B meson  mode  the  antineutron  shower 
energy  spectrum  is  dominantly  in  the  range  1.0  Emain  1*^  GgV. 
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Figure  3.9:  Emain  in  GeV  without  proton  requirement  for  n’s,  7’s  from  7t°’s,  and 
KEs 


Figure  3.10:  Ej^ain  in  GeV  with  proton  requirement  for  n s,  7 ^ from  tt  s,  and 
KEs 


CHAPTER  4 

MEASUREMENT  OF  ^ D*-~PP'k^ 

Modes  in  B [jD]AiVA  have  not  been  previously  measured.  The  reconstruc- 
tion efficiency  multiplied  by  the  number  of  B mesons  in  our  dataset  for  B^—^  D* 
p p and  D*~  p n , assuming  branching  fractions  similar  to  the  ones  mea- 
sured for  5+  ->  AcP7t+  and  ->  AcP7r+  7t“  [29],  is  at  a level  we  can  measure, 
which  lends  credibility  to  a search  for  both  modes. 

4.1  Monte  Carlo  Reliability 

The  reconstruction  process  at  CLEO,  as  well  as  at  all  other  high  energy  de- 
tectors, involves  a reconstruction  efficiency  toata-  For  s^riy  given  B decay,  only  a 
fraction  of  the  number  present  in  the  data  can  be  partially  or  fully  reconstructed. 
If  X number  of  B decays  in  our  data  are  B -4  Baryons  decays,  only  enata  x ^ 9-re 
measured.  In  the  case  of  R — > pX,  for  example,  there  is  a significant  contamina- 
tion due  to  protons  from  beam  gas  interactions  at  the  low  end  of  the  momentum 
spectrum  which  limits  B — > pX  to  be  carried  out  using  antiprotons  only,  and  then 
multiplying  the  result  by  2. 

By  tuning  a wide  range  of  parameters,  the  Monte  Carlo  simulation  can  be 
made  to  closely  resemble  the  data,  in  which  case  the  assumption  coata  ~ ^mc  Is 
reasonable,  cmc  is  reliable  because  several  processes  which  are  measured  at  CLEO 
allow  considerable  ease  of  simulation  as  well  as  large  samples  which  can  be  reliable 
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separated  from  backgrounds.  Processes  such  as  e'^e~,  e'^e~  — f 

e'^e~  — )■  7r“''7r~,  e'^e"  — >■  77,  and  e'^e“  — )•  777  allow  very  accurate  Monte  Carlo 
modeling  of  individual  particles,  which  in  turn  can  be  combined  to  model  many  B 
decays  well. 

On  this  level  of  accuracy  a second  group  of  decays  is  studied  to  widen  the  scope 
of  our  simulation:  Ks  TT'^Tr~,  (j)  — > D*  — >•  D^iisoft  with  , and 

A p7i  all  have  very  low  backgrounds.  Many  particle  separation  studies  are  based 
on  these  decays.  An  example  is  the  use  of  a low  background  sample  of  A’s  to  study 
antiproton  annihilation  showers,  as  will  be  discussed  later  in  this  chapter. 

4.2  Reconstruction  Procedure 

B mesons  are  produced  with  300  MeV  momentum  at  CLEO,  which  implies 
that  Ebeam  is  very  close  to  Eb-  This  momentum  is  small  when  compared  with  the 
B meson  mass  of  5.28  GeV.  Since  the  beam  energy  at  CLEO  is  measured  with 
a 2 MeV  resolution,  we  can  constraint  the  mass  of  the  B meson  candidate  to  be 
equal  to  the  beam  energy  Ebeam-  We  therefore  use  the  beam  contrained  mass  when 
reconstructing  R®— )•  D*~  p p : 

Mbc  = \/ Ebeam^  ~ Pi^ 

where  Ebeam  is  the  beam  energy,  on  average  5.29  GeV,  and  Xlf  Pi^  is  the  sum 
of  daughter  momenta  squared. 

The  energy  difference  between  the  beam  energy  and  the  energy  of  the  recon- 


structed B candidate,  defined  by: 
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— ^beam  ~ -^reconstructedB 


is  centered  at  0 and  has  a Gaussian  width  similar  to  the  Mgc  Gaussian  width, 
whereas  background  events  are  much  more  likely  to  form  a random  distribution  in 
A^;.  Selecting  B candidates  using  Ae  cuts  based  on  the  Monte  Carlo  is  useful  in 
separating  signal  from  background.  As  modelled  by  our  Monte  Carlo  for 
D*~  p p TT'^  is  centered  near  zero  and  has  a Gaussian  width  a < 100  MeV,  which 
allows  us  to  separate  the  signal  from  backgrounds  that  differ  by  a miss-measured 
extra  pion.  The  methodology  we  use  in  measuring  D*~  p p is  analogous 
to  that  used  in  previous  exclusive  reconstructions  at  CLEO  [45,  46]. 


4.3  Monte  Carlo  Study 

In  order  to  find  the  detector  efficiency  for  the  reconstruction  of  the  modes  we 
are  searching  for,  we  generate  Monte  Carlo  in  which  one  of  the  B mesons  in  the 
event  is  forced  to  decay  to  the  mode  we  are  reconstructing.  We  refer  to  this  sample 
as  signal  Monte  Carlo.  In  this  signal  Monte  Carlo  B mesons  decay  according  to 
phase  space. 

We  also  use  a sample  of  generic  Monte  Carlo  BB  events  on  which  we  run 
the  D*~  p P analysis  code.  The  generic  sample  is  composed  of  the  B 
decays  that  have  previously  been  measured  as  well  as  randomnly  generated  events 
assuming  inclusive  momentum  distributions.  There  are  no  B [D]NNX  modes 
in  the  generic  Monte  Carlo  sample.  We  use  the  generic  Monte  Carlo  sample  to 
model  our  backgrounds  due  to  modes  of  the  we  are  not  reconstructing. 
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Applying  the  particle  selection  criteria  outlined  in  Chapter  3 for  the  >■  D* 
p p decay  daughters,  we  fit  the  distribution  to  a single  Gaussian  for  signal 
and  a 1st  order  Chebyshev  polynomial  for  background.  We  show  in  Table  4.1  and 
Table  4.2  the  result  of  these  fits. 


Table  4.1:  AE  fit  results  for  D*  p p signal  MC  in  CLEO  II 


Decay  mode 

Central  value  (MeV) 

a (MeV) 

AE  cut  (MeV) 

D^^K+7^- 

0.6 

11.8 

ib  35 

-2.3 

16.8 

di  50 

1.1 

9.8 

± 29 

Table  4.2:  AE  fit  results  for  D*  p p signal  MC  in  CLEO  II. 5 


Decay  mode 

Central  value  (MeV) 

a (MeV) 

AE  cut  (MeV) 

-1.0 

9.8 

± 29 

-3.1 

13.4 

± 40 

D ® 7T  ~ 7T  7T  ~ 

-1.6 

6.1 

± 18 

We  define  mode/dataset  specific  AE  cuts  which  are  applied  to  Mbc-  Th®  Mbc 
distribution  is  then  fitted  to  a single  Gaussian  to  determine  each  tMCi  shown  in 
Table  4.3  and  Table  4.4. 

Table  4.3:  Mbc  distribution  fit  results  for  > D*~  p p tt"*"  signal  MC  in  CLEO 
II 


Decay  mode 

^MC  (%) 

cr  (MeV) 

8.69  ± 0.21 

2.63 

4.48  ± 0.15 

2.86 

D®  — > iL'*' 7T  “ 7T  7T  “ 

3.89  ± 0.14 

2.65 
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Table  4.4;  Mbc  distribution  fit  results  for  D*  p p signal  MC  in  CLEO 
II.5 


Decay  mode 

^MC  (%) 

a (MeV) 

4.32  ± 0.22 

2.55 

1.88  ± 0.14 

2.60 

— >•  7T  “ 7T  7T  ~ 

1.05  ± 0.08 

2.50 

The  TTsoft  in  D*  ->  D'^Ttgoft  for  the  CLEO  II. 5 dataset  is  the  most  sensitive 
particle  to  errors  in  our  Monte  Carlo  simulation.  The  replacement  of  the  PTL 
detector  by  the  SVX  detector  decreased  by  approximately  a factor  of  2 the  recon- 
struction efficiency  for  the  TTgoftS  used  here,  and  may  have  introduced  a systematic 
uncertainty. 


4.4  Results  in  Data 

In  Figure  4.1  and  Figure  4.2  we  show  the  A.E  vs  Mbc  distribution  for  )■ 
D*~  p p in,  respectively,  CLEO  II  and  CLEO  II. 5.  In  Figure  4.3  and  Figure 
4.4  we  show  the  Mbc  distribution  for  D*~  p p in,  respectively,  CLEO  II 
and  CLEO  II. 5,  after  the  mode  dependent  AE  cuts  have  been  applied.  The  mode- 
by-mode  Mbc  distributions  are  combined  for  each  of  Figure  4.3  and  Figure  4.4, 
in  which  (dark)  events  are  from  continuum.  The  continuum  background  to  )• 
D*~  p p TT"*"  is  statistically  insignificant.  If  this  had  not  been  the  case,  continuum 
subtraction  would  have  been  neccessary  to  subtract  this  background  from  the  Mbc 


distribution. 
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beam  constained  mass  (GeV)  in  CLEO  II 


5.300 


Figure  4.1:  AE  vs  Mbc  distribution  for  D*  p p in  CLEO  II  ON 

resonance  data. 
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Figure  4.2:  AE  vs  Mbc  distribution  for  D*  p p tc'^  in  CLEO  II. 5 ON 
resonance  data 
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Figure  4.3:  Mbc  (m  GeV)  for  D*  p p in  CLEO  II 
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beam  constained  mass  (GeV)  in  CLEO  11.5 


Figure  4.4;  Mbc  (in  GeV)  for  D*~  p p tt+  in  CLEO  II.5. 
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In  Table  4.5  we  attempt  various  fits  to  the  combined  CLEO  II  and  CLEO  II. 5 
Mbc  distribution  for  D*~  p p 'k'^  \ 

1.  Allowing  Mbc  and  ombc  float. 

2.  Fixing  Mbc  = 5.28  GeV  and  allowing  ombc  float. 

3.  Allowing  Mbc  to  float  and  fixing  umbc  according  to  expectations  from  signal 
MC. 

4.  Fixing  both  Mbc  and  umbc  prescribed  in  2.  and  3.  above. 

The  result  of  these  fits  are  shown  in  Table  4.5.  The  variation  in  yield  is  small. 
Table  4.5:  Results  in  the  CLEO  II/II.5  data  of  various  fits  to  Mbc  for  R®— ^ D*~ 

p p •K'^ 


fitting  options: 

fit  results 

floating  values 

Mbc  (in  MeV) 
a (in  MeV) 
Fit  yield 

5, 280.7iS:lg 
2.05iS:|| 

on  [X9-H6.55 

OU.OZ_5  92 

fixed  Mbc 

a (in  MeV) 
Fit  yield 

2.4o:^:^^ 

01  00+7.34 

Oi.OO„6  54 

fijfed  a 

Mbc  (in  MeV) 
C7  (in  MeV) 
Fit  yield 

5, 280.51^;^^ 
2.65 

32.441^:^1 

fixed  Mbc  and  a 

Fit  yield 

32.981^;^^ 

B{B^^  D*~  ppTT-^)  X 10“^ 

6.611:^ 

D*  p p ) is  calculated  using: 
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D*  p p 'n'^  ) = 


FittedYield 

(2  X X f^Mc)cLEOii'^(‘^  xB^B^x  ^mc)cleoii.5 


The  product  {BB  x ^mc)clboii  ~ product  {BB  x (■mc)cleoii  b- 
number  of  events  for  Mbc  > 5.275  GeV  is  21  for  CLEO  II,  and  15  for  CLEO  II. 5. 
Fitting  each  dataset  separately  would  yield  a measurable  B{B^-^  D*~  p p ). 

In  table  4.6  we  quote  the  number  of  background  events  in  the  range  5.2  GeV 
< Mbc  < 5.275  GeV,  and  the  number  of  signal  events  for  Mbc  > 5.275  GeV.  The 
mode  with  the  largest  background  is  D^—^K'^7r~7r^.  Signal  region  contains  « 
10  % background  events. 


Table  4.6:  J5°— >■  D*  p p : number  of  events  found  per  mode 


Mode 

Background 

region 

CLEO  II 

Signal 

region 

CLEO  II 

Background 

region 

CLEO  II.5 

Signal 

region 

CLEO  II.5 

d^-^k+b- 

7 

6 

0 

6 

D^^K'^‘K~7r^ 

14 

10 

15 

6 

4 

5 

3 

3 
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Figure  4.5:  Mbc  for  D*“  p p in  data 
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We  find  the  A E distribution  for  D*~  p p tt'*'  , as  shown  in  Figure  4.6  to 
be  Gaussian,  as  expected  from  signal  MC.  We  also  find  the  momentum  spectrum 
for  -Kg  from  the  D*  decay,  as  shown  in  Figure  4.7,  to  agree  with  expectations  from 
signal  MC  and  to  be  very  soft.  In  Figure  4.7  the  solid  distribution  is  CLEO  II 
Monte  Carlo  and  the  dashed  distribution  is  CLEO  II/II.5  data.  Only  events  with 
Mbc  > 5.27  are  plotted  for  both  figures. 


MINUITx^FittoPlot  541  &11 

Aa.  dele  axis 

File:  Generated  Internally 

Plot  Area  Total/Fit  41.000  / 41.000 

Func  Area  Total/Fit  40.999  / 40.999 

Y^=  1.6  for  14-  Bd.O.f., 

Errors  Paratxilic 

Function  1:  Chebyshev  Polynomial  of  Order  1 
NORM  280.17  + 88.06 

CHEB01  0.36391  ± 0.3856 

Function  2:  Gaussian  (sigma) 

AREA  21.387  ± 6.384 

MEAN  -2.50699E-03  ± 1.9945E-03 

SIGMA  6.57031  E-03  ± 1.6709E-03 


11 -APR-2000  14:46 
Fit  Status  0 
E.D.M.  1.00 

C.L.=  99.6% 

Minos 

O.OOOOE+00  -t-  O.OOOOE+00 

OPOOOE-rOO  + O.OOOOE-rOO 

O.OOOOE-100  + O.OOOOE-TOO 

O.OOOOE-rOO  -r-  O.OOOOE-rOO 

O.OOOOE-rOO  -r  O.OOOOE-rOO 


-0.040  -0.020  0.000  0.020 

Delta  E in  GeV 


Figure  4.6;  AF  (in  GeV)  for  D*~  p p 7t+  in  CLEO  II/II.5  data 
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Soft  pion  momentum  (GeV) 


Figure  4.7:  tTs  from  D*  momentum  for  jB°— ^ D*  p p ■k'^  in  data  and  Monte  Carlo 
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4.5  Resonant  Substructure 

In  order  to  determine  whether  we  are  using  the  correct  €mc  to  measure 
D*~  p p TT'^  )j  we  need  to  know  if  we  are  measuring  > D*~  p p or  a mode 
with  intermediate  particles  decaying  to  the  same  combination  of  decay  daughters. 
We  therefore  search  for  possible  contributions  to  the  resonant  substructure  of 
D*~  p p TT'^  . Any  measurable  resonant  substructure  would  have  to  be  subtracted 
from  the  inclusive  measurement  of  B{B^—^  D*~  p p 'k^  ),  and  the  €mc  would 
likewise  have  to  be  adjusted. 

4.5.1  Two-body  Decay  and  Possible  Strong  Resonances 

We  search  for  two  types  of  resonances: 

1.  A heavy  charmed  baryon  decaying  strongly  to  p + D*~ . No  significant 
peaking  is  observed. 

2.  A resonance  of  the  virtual  W decaying  to  pp-n'^ , B^-^  D*~  p p being  a 
two  body  decay.  No  significant  peaking  is  observed. 

We  study  the  Mp^p^T^+  spectrum  in  data  to  check  if  it  is  consistent  with  phase 
space  or  two-body  decay.  For  the  two-body  sample  we  generate  MC  with  a fictitious 
heavy  particle  that  decays  to  pp7[~^  for  ^ D*~  p p . It  has  a mass  of  2.6 
GeV  and  a width  of  200  MeV. 

No  conclusive  evidence  is  found  for  a measurable  contribution  to  B{B^-^  D*~ 
p p TT'^  ) from  these  possible  contributions.  The  6mc  for  B^^  D*~  p p varies 
slightly  depending  on  the  resonance  substructure.  We  allow  a 5%  systematic  un- 
certainty to  account  for  these  variations. 
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4.5.2  A Baryon  Contributions  in  the  Form  of  — > D*  and  5°  — 

D*-pA° 


Only  the  non-resonant  J5°— > D*~  p p are  shown  in  table  4.5.2.  The 

CMC  does  not  vary  by  more  than  5 % for  any  of  the  three  possible  assumptions: 
> D*~  p p D*~pA'^^,  and  D*~pA^.  Any  contribution  from 

the  latter  two  would  imply  a lumping  of  events  in  the  1.0  to  1.3  GeV  range  for 
Mp;r+  or  Mp^+.  The  % makeup  of  Mp^+  and  Mp^+  is  nearly  identical  in  Monte 
Carlo,  and  we  quote  only  the  former.  There  are  41  events  in  data  in  signal  region. 

Table  4.7:  % composition  of  generated  (phase  space  MC)  and  events  reconstructed 
(data)  of  Mp,r+  and  Mp,r+ 


Mass 
in  GeV 

Signal 
MC  (%) 

Data 

A/p7T+ 

% 

Data 

Mp^+ 

% 

1.0  - 1.1 

7.1  ± 2.1 

0 

0.0 

2 

4.9 

1.1  - 1.2 

11.2  ± 0.8 

3 

7.3 

4 

9.8 

1.2  - 1.3 

11.9  ± 0.7 

8 

19.5 

7 

17.1 

1.3  - 1.4 

11.4  ± 0.6 

8 

19.5 

9 

22.0 

1.4-  1.5 

10.2  i 0.6 

11 

26.8 

7 

17.1 

1.5  - 1.6 

9.2  ± 0.6 

3 

7.3 

■ 6 

14.6 

1.6  - 1.7 

11.3  ± 0.6 

3 

7.3 

3 

7.3 

1.7-  1.8 

6.4  ± 0.6 

3 

7.3 

3 

7.3 

We  place  an  upper  limit  of  < 2 events  for  Mp,r+,  and  < 3 events  for  Mp,r+ 
as  contributions  from  D*~pA'^'^,  and  B^  -4  D*~pA^  respectively.  These 

upper  limits  are  educated  guesses,  and  do  not  have  any  effect  on  the  systematic 
uncertainty  of  our  measurement  of  B{B^-¥  D*~  p p 'k^  ),  since  our  quoted  value 
includes  the  resonant  substructure.  Our  quoted  systematic  uncertainty  allows  for 
the  possibility  of  A baryon  contributions  to  the  resonant  substructure  of  J5°— > D*~ 


p p . 
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4.6  Backgrounds 

We  use  generic  Monte  Carlo,  which  is  discussed  in  Section  3.2,  to  study  possible 
backgrounds  to  )•  D*~  p p . Several  characteristics  of  >■  D*~  p p 
result  in  expected  low  backgrounds  as  well  as  an  accurate  Monte  Carlo  simulation: 

1.  J5°— D*~  p p TT'^  has  a proton  and  an  antiproton  as  decay  daughters,  the 
baryon-antibaryon  constraint  is  strictly  applied. 

2.  Z)*’s  have  low  BB  and  continuum  backgrounds. 

3.  The  high  combined  mass  of  the  decay  daughters  is  an  added  suppressant  of 
continuum  backgrounds. 

The  AE  cut  suppresses  modes  of  type  > D*~  p p 7r+  X,  where  X is  any 
additional  number  of  neutral  and  charged  pions,  but  it  does  not  suppress  their 
combinatoric  backgrounds.  Since,  as  previously  mentioned,  B — > [D]NNX  decays 
are  not  accounted  for  in  the  generic  MC  sample  we  study,  it  is  reasonable  to  assume 
that  if  75  % of  the  background  in  the  range  5.2  GeV  < Mbc  < 5.275  GeV  is  not 
accounted  for  by  the  generic  MC  sample,  the  most  likely  source  of  this  background 
is  combinatoric  background  to  B [D]NNX  decays.  The  small  amount  of  non- 
B — > [D]NNX  background  is  shown  as  a solid  filled  region  in  Figure  4.8,  in  which 
filled  overlay  on  data  is  combinatoric  non-B  [D]A^iVX  background  from  generic 
MC. 
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Figure  4.8;  Mbc  data  and  generic  Monte  Carlo  for  D*  p p (in  GeV) 
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4.7  Systematic  Uncertainties 


In  Table  4.8  we  show  the  systematic  uncertainties  we  consider  to  contribute 
significant  errors  to  our  measurement  of  B{B^^  D*~  p p ).  D*~  p p 

has  an  average  of  6.6  tracks. 


Table  4.8:  Estimate  of  systematic  uncertainties  (in  %)  for  > D*  p p 'k^ 


Source 

Uncertainty  (in  %) 

branching  fractions 

0.6 

D*  branching  fraction 

1.4 

D*  reconstruction 

5.0 

Monte  Carlo  statistics 

5.0 

# of  BB^s 

2.0 

tracking  (1%/track) 

6.6 

PRLEV  proton  ID 

8.0 

A signal  contribution 

5.0 

A background  contribution 

5.0 

Phase  space  versus  two  body 

5.0 

TOTAL 

15.5  % 

CHAPTER  5 

MEASUREMENT  OF  ->  D*~PN 


As  mentioned  in  Chapter  4,  we  have  selected  the  B — > [D]NNX  modes  which 
are  expected  to  have  the  lowest  backgrounds  as  well  as  the  highest  reconstruction 
efficiency.  When  reconstructing  D*~  p n , however,  we  encounter  a stumbling 
block  not  encountered  in  the  reconstruction  of  D*~  p p by  virtue  of  the 
former  mode  having  an  antineutron-or  a neutron  in  the  case  of  its  charge  conjugate- 
as  one  of  its  decay  daughters.  Our  inability  to  separate  neutron  showers  from  their 
backgrounds  reduces  the  reconstruction  efficiency  for  R®— >■  D*~  ' p n by  50%.  The 
antineutron  reconstruction  efficiency,  however,  is  sufficiently  large  (in  the  range  of 
30  to  40%).  We  do  not  have  an  accurate  measurement  of  the  antineutron  energy, 
yet  we  do  have  a well  measured  direction  of  its  shower,  which  we  use  to  reconstruct 
>■  D*~  p n . The  reconstruction  method  we  use  is  analogous  to  that  used  in  an 
unpublished  reconstruction  of  B Jj'^Ki  by  CLEO  [47].  In  Table  5.1  we  outline 
our  antineutron  selection  criteria  which  we  derived  in  Chapter  3. 

Table  5.1:  Antineutron  shower  selection  criteria. 

Track-to-shower  match  level  type  3 
E90E25  < cut  1 
|(cos(^)l  < 0.71,  or  good  barrel 
Emain  > 500  MeV 
Egroup  > 800  MeV 
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5.1  Reconstruction  Procedure 

Instead  of  Mbc,  we  define  m^o,  in  which  we  set  AE  = 0 and  assign  the  missing 
energy  to  the  antineutron  using  the  directional  cosines  of  its  shower.  Since  the 
measured  shower  energy  for  an  antineutron,  even  after  summing  the  energy  in  the 
list  of  nearby  showers,  fails  to  match  the  total  energy  of  the  antineutron,  we  only 
use  the  electromagnetic  shower  energy  as  part  of  our  selection  criteria. 

The  reconstruction  steps  for  this  mode  are: 

1.  En=Eheam-ED*+p  IS  assigned  to  the  antineutron.  This  energy  difference  is 
the  equivalent  of  setting  AE  = 0,  or  Ebeam  = EreconstmctedB- 

2.  3-momentum  magnitude  of  antineutron  candidate  pn=^E^  — m?,  with  rUn 
= 0.9395  GeV. 

3.  3-momentum  magnitude  of  antineutron  pn  times  x,y,z  directional  cosines  of 
shower  energy  are  assigned  to,  respectively,  x,y,z  components  of  Pn- 

4.  4-momentum  of  B candidate  pso  = Pn  + Pd*  + Pproton- 

5.  Mass  of  B candidate  meo  = J Elgam  ~ Pb°  • 

As  in  the  D*~  p p reconstruction,  we  use  Monte  Carlo  samples  in 

which  the  B meson  decays  according  to  phase  space.  The  resonant  substructure 
of  D’^—¥K'^7t~TT^  is  same  as  in  D*~  p p . 
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5.2  p n in  Monte  Carlo 

In  reconstructing  D*~  p n we  can  also  be  reconstructing  D*~ 

with  ->•  p h.  The  5°  ->  D*~  with  Df  p n contribution  is  part  of  our 

signal,  yet  it  is  not  a mode  within  the  resonant  substructure  of  B 0-4  D*-  ' p n . 
It  is  a separate  mode  of  the  B meson  with  B{B^  -4  Df  D*~)  = 9.6  ± 3.4  x 10“^ 
[5].  -4  p h,  a Feynman  diagram  of  which  is  shown  in  Figure  5.1.  Yet  another 

contribution  which  can  be  in  the  Mp+n  distribution  is  B^  -4  D*'^  D*~  with  Df 
p fi,  which  is  a background.  B{B^  -4  D*~^  D*~)  = 2.0±0.7  x 10“^  [5].  The  effect 
of  either  on  5° -4  D*~  p n cannot  be  assesed  accurately  since  we  do  not  know 
B{D^  — > p fi).  If  we  assume  B{Df  -4  p n)  is  ~ 1%,  the  effect  on  5°-4  D*~  p n 
can  be  significant  enough  to  affect  our  measurement  of  B{B^-^  D*~  p n )■ 


C U 


Figure  5.1:  A Feynman  diagram  for  —>■  p n 
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The  Mp^n  distribution  should  not  peak  appreciably  in  D*~  p n , and  a 

narrow  Gaussian  peak  with  5 MeV  a in  —y  Df  D*~  with  —^pn.  The  Df 

p n Monte  Carlo  peak  for  B'^  — > Df  D*~  is  shown  in  Figure  5.2. 


Ds  Mass  (GeV) 

Figure  5.2:  Mp+n  in  GeV.  B^  — >•  Df  D*~  with  Df  —ypn 


In  Figure  5.3  we  show  the  smeared  Mp^n  distribution  resulting  from  recon- 
structing B^—^  D*~  p n in  a signal  Monte  Carlo  B'^  — >■  D*~^  D*~  with  -4  p 
n sample.  ->■  Dfj  = 1.  The  missing  soft  photon  will  cause  this  background 
sample  to  be  considerably  broader  than  the  B^  — > D*~  with  Df  p n signal 

contribution. 

The  D*~  contribution  is  a background  to  D*~  p n ns  well 

as  B^  Df  D*~  with  p n.  Our  choice  regarding  how  to  deal  with  these 

contributions  is  to  exclude  both  possible  contributions  {B^  — > D*~  with 
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Figure  5.3;  Mp+n  (in  GeV)  from  a reconstruction  of  D*  p n in  a signal 
Monte  Carlo  — >■  D*~  with  Df  p n 

p n and  D*'^  D*~  with  p h)  in  quoting  our  final  value  for  B{B^-^ 


D*  p n ). 


Both  Df  p n contributions  can  be  excluded  with  1.91  GeV  < Mp+n  < 2.04 


GeV.  In  Tables  5.2  and  5.3  we  outline  the  cj^c’s  for  D*  p n . These  cmc’s 


take  into  consideration  that  we  are  not  reconstructing  the  conjugate  mode-with  a 


neutron. 
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Ds  Mass  (GeV) 


Figure  5.4:  (white)  is  —¥  Df  D*  with  —^pn.  (dashed)  is  -4  D* 
with  Df  p n.  Mpj^n  (in  GeV) 


Figure  5.5:  (white)  is  — » Df  -D*  with  Df  p n.  (dashed)  is  B^  D*"*"  D* 

with  p n.  rriBO  (in  GeV) 
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Table  5.2:  )•  D*  p n signal  MC  e^c’s  and  widths  in  CLEO  II. 


Mode 

% 

MeV 

both  Df  p n contributions  excluded. 

D^-^K+7^- 

7.47 

ih 

0.28 

3.10 

2.88 

± 

0.17 

3.40 

D^-^K'^TT~ir'^TT~ 

3.33 

± 

0.19 

3.11 

both  Df  p n contributions  included, 

without  any  background  correction. 

8.12 

± 

0.32 

3.14 

D^-^K+TT-7T^ 

3.09 

± 

0.19 

3.00 

D'^—^K'^TT~n'^TT~ 

3.52 

± 

0.20 

3.05 

Table  5.3:  jB°-4  D*~  p n signal  MC  eMc’s  and  widths  in  CLEO  II. 5 


Mode 

% 

MeV 

Both  Df  p n contributions  excluded. 

4.51  ± 0.17 

2.95 

1.60  ± 0.13 

3.23 

D^-^K'^n''T:'^Tr~ 

1.68  ± 0.14 

2.79 

Both  Df  —^pn  contributions  included, 

without  any  background  correction. 

D^^K+tt- 

4.99  ± 0.18 

2.92" 

f)MJT+7T-7r° 

1.82  ± 0.14 

3.35 

D*^  — > K"*"  7T  “ 7T  7T  " 

1.83  ± 0.15 

2.78 
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5.3  Results  in  Data 

We  look  for  D*~  p n with  or  without  the  mass  region  1.91  GeV  < Mp+n 
< 2.04  GeV.  We  compare  in  Figure  5.6  these  two  cases  in  CLEG  II/II.5. 

16 
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5.200  5.225  5.250  5.275  5.300 

B mass  (GeV) 

Figure  5.6:  (white)  is  inclusive  of  all  contributions,  (dashed)  is  after  exclusion  of 

both  p n contributions 

« 

The  D*~  p n m^o  distribution  is  less  statistically  significant  without  the 

Df  —¥  p n contributions.  However,  the  cost  of  lost  events  for  the  purposes  of  our 

( 

measurement  is  warranted  for  the  following  reasons: 

1.  The  signal  events  we  are  left  with  have  an  insignificant  probability  of  being 

something  other  than  p n . 

2.  We  properly  account  for  the  lost  events  in  our 

In  Figure  5.7  we  show  m^o  for  ON  resonance  and  OFF  resonance  for  the  case 
we  use  to  quote  D*~  p n ):  without  both  contributions. 
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Figure  5.7:  (white)  is  ON  resonance,  (solid)  is  OFF  resonance 


The  number  of  OFF  resonance  events  is  not  significant  enough  to  affect  our 
measurement.  We  choose  not  to  subtract  these  from  the  ON  resonance  distribution. 
In  Figure  5.8  and  Figure  5.9  we  show  m^o  distributions  by  dataset. 
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Figure  5.8:  for  D*  p n in  CLEO  II  ON  resonance  data 
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Figure  5.9:  Mbo  for  D*~  pn  in  CLEO  II.5  ON  resonance  data 
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5.4  D+  —^pnin  Data 

In  Figure  5.10  we  plot  Mp+n  for  the  region  where  we  expect  to  see  Df  p n 
as  a small  width  signal  for  — > Df  D*~ , and  as  a shallow  background  for 

— > Df^  D*~ . The  horizontal  lines  demarcate  the  mass  region  1.91  GeV  < Mp^n 
< 2.04  GeV,  which  we  are  excluding.  Events  with  m^o  > 5.275  GeV  region  are 
shown  in  Figure  5.10.  The  vertical  lines  demarcate  the  mass  region  1.91  GeV  < 

A 

Mp+n  < 2.04  GeV,  which  we  are  excluding.  Not  shown  are  events  with  > 

2.27  GeV.  We  expect  to  find  10%  of  B^-^  D*~  p n signal  events  in  this  region, 
which  translates  to  3 events. 


mass  of  proton  + antineutron  in  CLEO  II/II.5 


Figure  5.10:  Mp+n  for  >•  D*  p n (in  GeV)  in  CLEO  II/II.5 

We  find  8 events  where  a Df  p n signal  is  expected.  This  sample  is  not 
sufficiently  significant  to  label  these  events  as  signal.  The  scatter  plot  in  Figure 
5.11  contains  more  information  than  Figure  5.10:  there  are  other  regions  in  this 
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plot  that  are  equally  statistically  significant  to  the  region  where  we  expect  to  find 
Dt^  p n.  If  we  are  to  claim  a signal  in  one,  we  should  be  able  to  claim  one  in 
any  of  these  others. 

2.67 
2.47 
2.27 
2.07 


1.87 

5.200  5.225  5.250  5.275  ' 5.300 

B mass  (in  GeV)  in  CLEO  II/II.5 

Figure  5.11:  Mp+s  vs.  m^o  (both  in  GeV) 

We  can,  however,  estimate  an  upper  limit  for  D+  ph.  The  most  conservative 
estimate  is  to  call  all  8 events  in  the  expected  region  signal,  which  leads  to  a value 
of  11.5  for  a 90%  confidence  level  measurement.  11.5  signal  events  correspond  to 
->  p h)  < 5.4  % using  the  PDG  value  for  D*~). 
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5.5  D*  p n ) Measurements. 

I 

We  use  the  same  equation  to  calculate  B{B^^  D*~  p h ) as  we  did  for  > 
D*~  p p TT'^  : 


B{B^-^  D*  p p ) = 


FittedYield 

(2  X X ^Mc)cLEOii~^('^  ^ B^B^  x ^mc)cleoii  .h 


As  we  mentioned  when  calculating  B{B^^  D*  pp  ),  the  product  {BB  x cmc) 
is  almost  the  same  for  CLEO  II  and  CLEO  II. 5.  The  number  of  events  for  this 
mode  for  Mbc  > 5.275  GeV  is  10  for  CLEO  II,  and  16  for  CLEO  II.5. 

The  combined  m^o  distribution  can  be  fitted  to  yield  results  that  are  consistent 
with  Monte  Carlo  expectations. 


B mass  (GeV) 
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Figure  5.12:  Inclusive  m^o  in  data  for  > D*  p n (in  GeV) 
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Figure  5.13:  m^o  in  data  for  D*  p n in  GeV  excluding  both  p n 

contributions 


Table  5.4:  Measurements  of  D*  p n branching  fractions 


Mode 

Yield 

. Branching  Fraction 

B^^  D*~  p n inclusive 
B^-^  D*~  p n without  both  contributions 

32.941^-^^ 

24.02lti° 

14.9l2:7  X 10"^ 
11.9111  X 10~^ 

The  difference  in  branching  fractions  does  not  imply  that  what  has  been  ex- 
cluded is  only  ->  p h.  We  do  not  attempt  to  derive  an  upper  limit  by  this 
method.  As  in  the  case  of  our  measurement  of  )■  D*~  p p , we  find  the 
D*~  p n signal  events  are  consistent  with  phase  space  decay.  We  look  for  a 
resonance  of  type  (cud)  -I-  {dd)^^^  h -f  and  fail  to  find  it. 

The  fits  are  quoted  in  Table  5.5.  m^o  is  fitted  to  a single  Gaussian  for  signal, 
and  an  Argus  function  for  background.  The  offset  is  fixed  at  0.0025,  and  the  Ebeam 
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is  fixed  at  5.29  GeV  in  all  the  fits.  As  in  the  case  of  >■  D*~  p p , in  table 
5.5,  leaving  some  parameters  fixed,  but  not  others  allows  us  to  test  the  validity  of 
our  results.  The  excluded  region  is  1.91  GeV  < Mp+a  < 2.04  GeV. 

Table  5.5:  Results  of  various  fits  to  meo  for  D*~  p n with  or  without  the 
Df  —ypn  contributions 


Fitting  options 

No  p n 

Inclusive 

floating  values 

rriBO  (MeV) 
o (MeV) 
Fit  Yield 

5, 280.31^1 
2.59lg;^? 

OQ  p;7+5.61 
Zo.O/_4  97 

5,280.5i^i^ 

9 /1 9+0-46 

SlASttil 

fixed  ruBO  — 5.28  GeV 

a (MeV) 
Fit  Yield 

9 ^9+u.4y 
Z.OZ_o.42 

23.671^;^ 

9 p;p;+0-47 
Z.00_o.40 

31.97t®|^ 

fixed  a = 3.07  MeV 

rrigo  (in  MeV) 
Fit  Yield 

5, 280.31^:^7 
23.971^96 

5, 280.31^:^6^ 

99  7C7+6.56 

OZ.  / 0_5  91 

fixed  m^o  and  a 

Fit  Yield 

24.02+t^^ 

32.94i^;^^ 

S(R°-^  D*~  pn)  X 10-^: 
From  raw  yield 

11.9lli 

14.911;? 

With  MG  Correction  Factor 

u.shli 

ISAtli 

In  the  following  section  we  discuss  the  correction  factor  we  apply  to  the  raw 
D*~  p n)  due  to  a discrepancy  between  data  and  Monte  Carlo  in  the  ep  for 
antiproton  annihilation  showers.  The  Correction  Factor  we  use  increases  the  raw 
D*~  p n ) quoted  in  Table  5.5  by  a factor  of  1.21  with  an  added  systematic 
uncertainty  of  4%. 
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5.6  Correction  Factor 

The  Monte  Carlo  we  use  has  not  been  optimized  to  model  nucleon-antinucleon 
annihilation.  We  find  a discrepancy  between  the  reconstruction  efficiency  for  data 
and  Monte  Carlo  for  antiprotons.  We  consider  the  discrepancy  credible  and  change 
the  antineutron  efficiency  using  a Correction  Factor.  We  assume  that  the  Monte 
Carlo  fails  to  model  antineutrons  by  the  same  amount  as  it  does  for  antiprotons, 
even  though  the  quark  content  of  antiprotons  and  antineutrons  is  different. 

We  define  CannihUation  ^s  the  efficiency  for  an  antibaryon  shower  to  pass  the 
annihilation  shower  cuts  outlined  in  Table  3.8.  This  quantity  will  be  different  for 
antiprotons  and  antineutrons.  In  the  former  case  we  choose  a sample  which  is 
inevitably  biased  by  the  momentum  range  in  which  we  can  separate  antiprotons 
from  other  charged  tracks.  In  the  latter  case  we  use  a sample  which  has  shower 
backgrounds. 

We  define  a correction  factor  to  account  for  the  discrepancy  between  data  and 
Monte  Carlo  for  annihilation  showers  as; 


Correction  Factor  = 


^annihilation  ^MC 
^annihilation,  Data 


This  Correction  Factor  is  found  for  antiproton  momentum  bins  in  a range 
consistent  with  the  momentum  range  of  antineutrons  in  D*~  p h as  found 
from  the  energy  assigned  to  the  antineutron  candidate  shower  after  assuming  the 
antineutron  mass.  The  number  of  antineutrons  in  each  of  these  momentum  bins 
is  used  to  weigh  the  contribution  of  each  bin  to  the  Correction  Factor,  which  is 
defined  for  the  entire  range. 


no 


5.7  Use  of  a A Sample 

The  selection  of  a clean  sample  of  antiprotons  is  achieved  by  reconstructing  A’s. 
In  the  decay  A — > due  to  the  large  difference  in  mass  between  the  p and  the 

TT'*',  the  kinematics  of  the  decay  allow  for  a clear  separation  of  the  p tracks  from 
the  TT'^  tracks.  Furthermore,  the  single  Gaussian  width  of  the  distribution  is  1 
MeV  and  is  affected  by  backgrounds  insignificantly. 

Antiprotons  are  selected  with  the  same  criteria  used  in  the  reconstruction  of 
J5°-^  D*~  p p except  for  the  use  of  Lproton  > 0.1.  A’s  are  from  the  KNVF 
package  [42].  The  A selection  criteria  is  as  shown  in  Table  5.6. 

Table  5.6:  A selection  criteria 

Both  tracks  pass  TRKMNG  flag  tng(track)  > 0 
of  two  tracks  to  form  a vertex  < 30 
for  A to  point  back  to  the  interaction  region  < 30 
Significance  of  the  two  dimensional  flight  distance  < 3a 

Flight  distance  > 0.005  meter 

We  count  the  number  of  A’s  in  the  signal  region  after  applying  a double  Gaus- 
sian fit  with  fixed  parameters  as  shown  in  Table  5.7. 

Table  5.7:  A double  Gaussian  fixed  parameters 


Mass 

1.1158  GeV 

a^ 

0.00162 

Area2  / areal 

0.46 

^MEAN 

0 

a^a^ 

0.425 

Ill 


In  Table  5.8  we  find  the  % of  antiproton  showers  in  a A sample  which  pass 
the  annihilation  shower  cuts.  In  Table  5.8  we  find  ^annihilation  for  MC  and  Data  A 
samples.  The  samples  are  a combination  of  ON  and  OFF  resonance.  The  ON  and 
OFF  resonance  e annihilation^  were  found  to  yield  the  same  results. 


Table  5.8:  e annihilation-  (%)  of  match  1 or  2 antiprotons  in  A’s  passing  annihilation 
shower  selection  cuts  for  six  momentum  spectra 


Momentum  range 

CLEO  II 
MC 
% 

CLEO  II 
Data 
% 

CLEO  II.5 
MC 
% 

CLEO  II.5 
Data 
% 

300-500  MeV 

28.3  ± 1.3 

13.2  ± 0.3 

26.5  ± 1.0 

13.3  ± 0.2 

500-700  MeV 

65.1  ± 1.7 

52.7  ± 0.9 

65.0  ± 2.6 

52.2  ± 0.6 

700-900  MeV 

70.8  ± 1.8 

60.6  ± 1.0 

71.5  ± 2.8 

60.9  ± 0.7 

900-1100  MeV 

65.3  ± 1.7 

52.9  ± 0.9 

63.4  ± 2.5 

50.7  ± 0.6 

1100-1300  MeV 

50.6  ± 2.0 

41.9  ± 1.0 

52.2  ± 2.1 

41.4  ± 0.7 

1300-1700  MeV 

47.4  ± 1.9 

44.0  ± 1.0 

46.3  ± 1.8 

41.0  ± 0.7 

1700-2100  MeV 

56.0  ± 2.0 

54.0  ± 1.0 

55.1  ± 1.9 

53.2  ± 0.7 

The  Correction  Factor  is  weighed  according  to  the  fraction  of  antineutrons  in 
a momentum  range  as  produced  in  Monte  Carlo.  The  solid  distribution  in  Figure 
5.14  is  data  and  the  dashed  distribution  is  signal  Monte  Carlo. 
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Figure  5.14:  Antineutron  momentum  distribution  (in  GeV) 
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In  Table  5.9  we  calculate  the  Correction  Factor  as  previously  defined,  as  well 
as  the  fraction  of  antineutrons  found  in  signal  Monte  Carlo  for  each  momentum 
range. 


Table  5.9:  Correction  Factor  (C.F.)  as  a function  of  Pprotm  ^-nd  % of  antineutrons 
found  in  signal  Monte  Carlo  for  each  momentum  range  in  CLEO  II/II.5 


Momentum  range 

% of  n in  MC 

C.F. 

300-500  MeV 

9.1 

2.07  ± 0.11 

500-700  MeV 

10.0 

1.25  ± 0.04 

700-900  MeV 

14.7 

1.17  ± 0.04 

900-1100  MeV 

15.2 

1.24  ± 0.04 

1100-1300  MeV 

16.3 

1.23  ± 0.06 

1300-1700  MeV 

20.7 

1.10  ± 0.05 

1700-2100  MeV 

15.7 

1.04  ± 0.05 

300-1700  MeV 

100.0 

1.21  ± 0.05 

We  will  use  1.21  as  the  Correction  Factor  by  which  to  increase  D*~  p 

n ),  which  is  an  equally  weighed  average  of  CLEO  II  and  CLEO  II. 5,  and  add  to 
our  list  of  systematic  errors  a 4%  contribution  due  to  this  correction. 
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5.7.1  Backgrounds  in  Signal  and  Generic  Monte  Carlo 

To  check  that  our  antineutron  selection  criteria  does  not  contaminate  our  sam- 
ple with  showers  from  decays  similar  to  > D*~  p n , we  run  our  D*~  p n , 
with  reconstruction  code  on  signal  Monte  Carlo  for  > D**  p n with 

equal  numbers  of  Di(2420)“  and  D2*(2460)°,  D*~pn7:^,  and  B°  -4  D*~pA^ 

with  A°  -4  n7T°.  The  result  for  m^o  are  overlayed  in  solid  on  our  data  results,  in 
white,  for  B°-4  D*~  pn  in  Figure  5.15. 


300 


200 


100 
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5.200  5.225  5.250  5.275  5.300 

Figure  5.15:  m^o  in  GeV  for  D*~  p n selection  criteria  applied  to  selected 
signal  Monte  Carlo  background  modes 

We  also  run  the  >•  D*~  p n code  for  all  three  modes  on  our  generic 
Monte  Carlo  sample,  as  shown  in  Figure  5.16.  The  generic  Monte  Carlo  sample  is 
discussed  in  Section  3.2.  Scaling  the  generic  Monte  Carlo  sample  to  the  size  of  our 
dataset,  we  find  it  to  account  for  40%  of  the  background. 


r 1- 


T 1 1 r 


1“ 


115 


Within  the  generic  MC  sample  ~ 25%  is  combinatoric  background  events  from 
each  of  D^-^K'^'K~T{'^Tr~  and  . The  remaining  events,  50%,  are  from 

D^^K'^n~7T^.  Since  B D*^NNX  decays  are  not  included  in  the  generic  MC, 
these  events  are  from  random  combinations  of  a 5 — > Baryons  decay  and  & B 
D*  X decay. 


10 
0 

5,200  5.225  5.250  5,275  5.500 

Figure  5.16:  rngo  for  B^-^  D*~  p n in  GeV.  (white)  is  data,  (solid)  is  generic  MC 

We  are  unable  to  account  for  60%  of  the  background  found  in  data  in  the  m^o 
distribution  for  B^^  D*~  p n . However,  this  background  is  small,  and  the  m^o 
distribution  we  find  is  consistent  with  a signal  at  the  nominal  B mass. 

5.8  Antineutron  Directional  Cosine  Resolution 

We  compare  the  generated  and  reconstructed  directional  cosines  for  the  an- 
tineutron candidates  in  signal  Monte  Carlo.  We  find  the  median  error  to  be  37.8 
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milliradians.  Since  a single  5-cm  calorimeter  cell  corresponds  to  w 50  milliradians, 
the  systematic  error  in  our  measurement  of  the  directional  cosines  of  the  antineu- 
tron candidate  shower  is  insignificant. 

5.9  — )■  D*^pn. 

Neutrons  do  not  annihilate  in  the  calorimeter.  Since  our  selection  criteria  for 
antineutron  showers  rejects  the  vast  majority  of  neutrons,  we  derive  a new  set  of 
cuts  specific  to  neutrons.  The  selection  criteria  we  used  is  shown  in  Figure  5.10. 
The  most  important  diflference  with  antineutrons  is  the  selection  of  low  energy 
photon-like  showers. 


Table  5.10:  Neutron  shower  selection  criteria 


Track-to-shower  match  level  type  3 
E90E25  > cut  1 
|(cos(0)|  < 0.71,  or  good  barrel 


E. 


main 


< 500  MeV 


We  are  unable  to  define  selection  criteria  that  allow  us  to  separate  neutrons 
from  soft  photons.  This  background  overwhelms  the  signal,  making  unatainable 
the  reconstruction  of  — )■  D*'^pn.  We  plot  Emain  for  neutrons  (solid  distribution) 

and  antineutrons  (dashed  distribution)  in  Figure  5.17. 
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Figure  5.17;  Emain  (in  GeV)  for  neutrons  and  antineutrons  in  D*  pn  signal 
Monte  Carlo 
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5.10  Systematic  Uncertainties 


In  Table  5.11  we  show  the  systematic  uncertainties  we  consider  to  contribute 
significant  errors  to  our  measurement  of  D*~  p n).  D*~  p n has  an 

average  of  4.6  tracks. 


Table  5.11:  Estimate  of  systematic  uncertainties  (in  %)  for  B^-^  D*  p n 


Source 

Uncertainty  (in  %) 

branching  fractions 
D*  branching  fraction 
D*  reconstruction 
Monte  Carlo  statistics 
# of  BB's 
tracking  (1  %/track) 
PRLEV  proton  ID 
A background  contribution 
Phase  space  versus  two  body 
antineutron  miss-ID 
TOTAL 

0.6 

1.4 

5.0 

5.0 

2.0 
4.6 

4.0 

5.0 

3.0 

15.0 
18.8  % 

CHAPTER  6 
CONCLUSION 


This  work  comprises  the  second  successful  exclusive  reconstruction  of  B ^ 
Baryons  modes,  after  the  exclusive  reconstructions  of  AcPTr"^  and  B^ 

AcPTT'*'  7T“  carried  out  by  CLEO  in  1997  [29].  In  Table  6.1  we  summarize  all  the 
exclusive  B — >•  Baryons  measurements  to  date.  We  rank  the  modes  from  largest 
to  smallest  branching  fraction  central  value.  The  first  error  is  statistical  and  the 
second  is  systematic.  The  Ac  modes  have  a third  uncertainty,  which  is  systematic 
and  is  due  to  the  error  in  the  determination  of  B(A^  — > pK~ir'^). 

« 

Table  6.1:  Exclusive  measurements  in  B — >■  Baryons  to  date 


Mode 

Branching  Fraction  (xlO  ‘^) 

B^-^  D*~  p n 
B^  — >•  KcjyK'^  7T“ 

B^^  D*~  p p TT"*" 

B~^  AcPTr"*" 

14.5ty  ± 2.7 
13.3lt|  ± 3.1  ± 2.1 
6.6^1;^  ± 1.0 
6.2l|^±l.l±1.0 

As  shown  in  Table  6.1,  the  size  of  the  branching  fractions  we  are  measuring 
in  this  work  are  of  the  same  order  of  magnitude  as  those  previously  measured  in 
decay  modes  including  a Ag.  We  are  unable  to  extrapolate  the  exclusive  results 
to  compare  the  relative  magnitude  of  the  inclusive  modes  with  or  without  a Ac. 
For  instance,  the  available  phase  space  for  the  production  of  more  light  mesons 
is  larger  in  the  case  of  modes  with  a Ac  than  in  the  case  without,  which  may 
cause  the  B — >■  AcX  inclusive  rate  to  be  larger  than  the  B — )■  [D]NNX  inclusive 
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rate.  However,  we  think  it  reasonable  to  assume  that  the  measurements  we  have 
at  our  disposal  to  date  are  representative  of  the  inclusive  modes.  We  infer  that 
B — > [D]NNX  modes  contribute  significantly  to  the  total  B — )■  Baryons  decay 
rate. 


6.1  B Baryons  phenomenology 

The  B Baryons  theoretical  models  attempted  to  date  [48,  49]  rely  on  the 
assumption  that  the  B meson  decays  to  two  decay  daughters,  which  is  the  two- 
body  assumption.  Unlike  two-body  B meson  modes,  most  B Baryons  modes, 
including  all  four  that  have  been  measured  to  date,  have  at  least  three  decay 
daughters.  Most  B Baryons  modes  also  have  a varied  resonant  substructure. 
The  assumption  that  the  hadronization  process  takes  place  at  a late  stage  in  the 
decay,  which  allows  for  a substantial  simplification  of  the  equations  due  to  the 
suppressed  dependence  on  the  exchange  of  gluons  and  light  quarks,  is  justified  in 
hard  two-body  decays.  However,  in  many-body  decays  there  can  be  re-scattering 
of  the  hadronizing  quarks  and  multiple  exchanges  of  soft  gluons  and  virtual  quarks. 
Notwithstanding  the  calculational  difficulties  of  many-body  B — ^ Baryons  decays 
such  as  the  ones  we  measure  here,  progress  continues  to  be  made  using  HQET  to 
explain  B meson  decay.  The  focus  has  been  on  decays  for  which  HQET  can  be  used 
to  extract  the  most  information,  and  which  can  be  used  to  search  for  CP  violation 
[14,  15].  Future  attempts  to  theoretically  explain  B Baryons  are  needed. 
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6.2  Possible  Future  B Baryons  Modes  at  CLEO 

In  CLEO  II  reconstruction  from  B decays  is  plagued  by  high  backgrounds 
due  to  poor  separation  between  kaons  and  pions.  The  substantial  improvements  in 
charged  particle  separation  (pion-kaon-proton)  in  the  CLEO  III  data  should  allow 
for  the  successful  reconstruction  of  B'^  — > D^pn  and  B~^  D^ppir'^. 

The  ability  to  reconstruct  antineutron  modes  should  allow  for  the  reconstruc- 
tion of  modes  of  type  B (Ac,Ec)nX.  Modes  of  type  B -4-  (Ac* ,T,c*)pX , and 
B HcAc,  which  have  low  reconstruction  efficiencies  in  CLEO  II/II.5,  are  also 
worth  pursuing. 

B —4  Baryons  for  b sg,  and  b ^ u modes  will  also  be  of  interest  and  likely 
within  range.  If  previous  measurements  of  the  meson  modes  are  used  as  a guide, 
the  modes  B^  -4  Ap,  B^  —4  pnTr~ , and  B^  —4  pp7r'^7r~  are  worth  pursuing. 

Semileptonic  b ^ c B Baryons  decays,  such  as  B^  -4  Ac~ pl'^i^i,  with  I — e,  p, 
due  to  the  low  reconstruction  efficiency  at  CLEO  for  low  momentum  leptons,  have 
not  been  measured.  The  substantially  larger  dataset  expected  in  CLEO  III  may 
compensate  for  the  low  reconstruction  efficiency.  A mode  such  as  B^  -4  ppl'^vi 
should  be  feasible.  This  latter  mode  can  be  used  as  an  auxiliary  in  the  measurement 

of  lKb|. 

6.3  Significance  of  Results 

We  have  found  the  first  evidence  of  decay  modes  of  the  B^  of  the  type  B — 4 
DNNn.  We  measure  the  branching  fractions  B (R°-4  D*~  p p tt'*')  = (6.6li;2  ± 
1.0)  X 10“^,  and  B{  R°-4  D*~  pn)  = (14.8t|;i  ± 3.0)  x 10“^.  These  measurements 


indicate  the  fraction  of  baryonic  decays  of  B mesons  that  do  not  proceed 
may  be  of  approximately  the  same  magnitude  as  those  that  do. 
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